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Abstract 

Waves in extended periodic structures are well-known to spatially disperse and decay in am- 
plitude as time advances. This dispersion is associated with the continuous spectrum of the 
underlying differential operator and the absence of discrete eigenvalues. The introduction of lo- 
calized perturbations, leads to defect modes, states in which energy remains trapped and spatially 
localized. In this paper we present results on the bifurcation of such defect modes, associated with 
the emergence of discrete eigenvalues from the continuous spectrum. 

1 Introduction 

Waves in extended periodic structures are well-known to spatially disperse and decay in amplitude 
as time advances. This dispersion (Floquet-Bloch dispersion) is associated with the continuous 
spectrum (extended states) of the underlying differential operator and the absence of discrete 
eigenvalues (localized bound states). The introduction of localized perturbations leads to defect 
modes, states in which energy remains trapped and spatially localized. This phenomenon is of great 
importance in fundamental and applied science - from the existence of stable states of matter in 
atomic systems to the engineering of materials with desirable energy transport properties through 
localized doping of ordered materials. 

The process by which the system undergoes a transition from one with only propagating delo- 
calized states to one which supports both localized and propagating states is associated with the 
emergence or bifurcation of discrete eigenvalues from the continuous spectrum associated with the 
unperturbed periodic structure. In this paper, we discuss this bifurcation phenomenon in detail 
for the Schrodinger operator 

Hq = -d'i + Q(x) (1.1) 
where Q(x) is a real- valued, sufficiently smooth periodic potential 

Q(x + 1) = Q{x). (1.2) 

The spectrum, spec(ifQ), of the Schrodinger operator is continuous and is the union of closed 
intervals called spectral bands |25j . The complement of the spectrum is a union of open intervals 
called spectral gaps. The spectrum is determined by the family of self-adjoint eigenvalue problems 
parameterized by the quasi-momentum k 6 (—1/2, 1/2]: 



Hqu(x;k) = Eu(x;k) , 
u(x + 1; k) — e 2nlk u(x; k) . 



(1.3) 
(1.4) 



That is, we seek fc— pseudo-periodic solutions of the eigenvalue equation. For each fc G (—1/2, 1/2], 
the self-adjoint eigenvalue problem (1.3 1-(1.4 1 has discrete eigenvalue-spectrum (listed with mul- 
tiplicity): 

E {k) < Ei(k) < •■• < E b (k) < ... (1.5) 
with corresponding k— pseudo-periodic eigenfunctions: 

u b (x;k) = e 2n ' kx p b (x;k), p b (x + l;k) = p b (x;k), b > 0. (1.6) 



The b spectral band is given by 



The spectrum of Hq is given by: 



B b = [J E b (k). (1.7) 

fce(-i/2,i/2] 



spec(// Q ) = [J B 6 = |J |J Sj(fc). 

6>0 6>0 fce(-l/2,l/2] 



Since the boundary condition is invariant with respect to k i— >• fc + 1, the functions E b (k) 

can be extended to all R as periodic functions of fc. The minima and maxima of E b (k) occur at 
fc = fc* € {0, 1/2}; see Figure [l] At such values of fc*, 9j?-E(, is either strictly positive or strictly 
negative; see Lemma |2.2| 

Consider now the perturbed operator Hq+v, where V(x) is spatially localized. By Weyl's 
theorem on the stability of the essential spectrum, one has spec cont (-ffq+v) = spec cont (-ff<2) |25| . 
The effect of a localized perturbation is to possibly introduce discrete eigenvalues into the open 
spectral gaps. Note that Hq+v does not have discrete eigenvalues embedded in its continuous 
spectrum; see [261 114| . 

In this paper we present a detailed study of the bifurcation of localized bound states into gaps 
of the continuous spectrum induced by a small and localized perturbation of Hq: 

H Q+XV = -d 2 x + Q(x) + XV(x) . (1.9) 

Here Q(x) is a sufficiently smooth 1— periodic function defined on H and V(x) is spatially localized. 
The parameter A is assumed to be positive and will be taken to be sufficiently small. We next 
turn to a summary of our results. See Theorem |3 . 1 1 and Theorem |3.4| for detailed statements. 

Let E+ = E bt (k t ), fc„ £ {0, 1/2}, denote an endpoint (uppermost or lowermost) of the (6*)" 1 
spectral band, bordering a spectral gap. We show that under the condition: 

d 2 k E bt (k,) x f \u bt {x;k,)\ 2 V(x) dx < 0, (1.10) 
Jo 

the following holds: There exists a positive number, Ao, such that for all < A < Ao, Hq + w has 
a simple discrete eigenvalue 

E(X) — E t + X 2 fi + 0(X 2+a ), a>0. (1.11) 

which bifurcates from the edge, E+ = E bt (k*), of B bt into a spectral gap. 

1. If dlE bt (k,) > and \u b , (x; fc»)j 2 V(x) dx < 0, then fx < and E(X) lies near the 
lowermost edge of B bt ; see the center panel of Figure [l] 

2. If dlE bt (k,) < and \u b , (x; fc»)| 2 V(x) dx > 0, then fx > and E(X) lies near the 
uppermost edge of B bt \ see the right panel of Figure [l] 
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For < A < Ao, ip x (x), the eigenstate corresponding to the eigenvalue, E(X), is well-approximated 
in L°° by, go(\x), where go(y) denotes the unique eigenstate of the homogenized operator 



H, 



i>»,efl = 



d d 

■~T Ab,,eB ~y- + Bb,,eH 0(y) , 

dy dy 



with constant effective parameters A&^eff an d Bi>,,eflr- Here, 

1 



„ ^diEb,{k*) 



^6,, eft' = 

on - 

is the inverse effective mass associated the the spectral edge, E* — Eb, (k* 

6«, c ff = / \ub, (x; fc»)| 2 V(x) dx , 



Si 



(1.12) 



(1.13) 



(1.14) 



and S(y) denotes the Dirac delta mass at y = 0. 

Remark 1.1. In the case where, Q = 0, then Ho = — d 2 and its spectrum consists of a semi-infinite 
interval, spec(_ffo) = [0, oo), the union of touching bands with no finite length gaps. Furthermore, 
Pb(x; k) = 1, for all \k\ < 1/2 and b > 0. The only band-edge band edge is located at E+ = Eo(0) = 
0, where we have: fc* = 0, E (k) = 4n 2 k 2 and dlE (k,) = 8ty 2 . 

In this case, our results describe the bifurcation of an eigenvalue from the edge of the continuous 
spectrum of Ho induced by a small and localized perturbation: H\v = —d 2 -\-W, under the condition 
J m V < 0; see the discussion below of }27f . 
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Figure 1: Sketch of spectra. Eigenvalues, Eb(k), k £ (—1/2, 1/2], b = 0, 1, 2, are displayed in green. The continuous 
spectrum, is in blue, and discrete eigenvalues are indicated through cross markers. Left panel corresponds spcc(Hq), Q 
periodic. The center panel corresponds spec(HQ_i_\y), where XV is small, localized negative. The right panel corresponds 
to spec(HQ+\y), where XV is small, localized and positive. 



1.1 Previous related work 

The case Q = 0, where Hq + w = — A + V(x) was considered by Simon |27| in one and two spatial 
dimensions (in spatial dimensions n > 3, it is well known that for sufficiently small A, —A + XV 
does not have discrete spectrum for all A sufficiently small). In one dimension, it is proved that if 
V is sufficiently localized and — oo < J V < 0, then H\y has a small negative eigenvalue E(X) of 



order A ; see the Corollary 3.3 and the discussion following it. Borisov and Gadyl'shin [T] study the 



spectrum of perturbations of one-dimensional periodic Schrodinger operators and obtain results 
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closely related to the current work. Their proofs rely on ODE techniques. In contrast, we use 
frequency space arguments, which can be extended to higher dimensions. Deift & Hempel [5] 
obtained results on the existence and number of eigenstates in spectral gaps for operators of the 
general type H — XW , where H has a band spectrum and W is bounded. Figotin & Klein [121 
113] studied localized defect modes in context of acoustic and electromagnetic waves. A formal 
asymptotic study, in terms of a Floquet-Bloch decomposition, in one and two spatial dimensions 
was given in Wang et al. |28j . Johnson et. al. |23) formulate a variational principle for defect 
modes with frequencies in spectral gaps. They use formal trial function arguments to show the 
existence of such defect modes in spatial dimensions one and two. By formal asymptotic arguments, 
they obtain the condition (1.10 1, for the case of the first spectral gap. Results on bound states 
and scattering resonances of one-dimensional Schrodinger operators with compactly supported 
potentials appear in work of Bronski & Rapti [5] and Korotyaev [191 .20] . Bifurcations of defect 
modes into spectral gaps was considered in dimensions d — 1, 2 and 3 by Hoefer & Weinstein |15j 
for operators of the form —A + Q(x) + e 3 V(ex), where Q is periodic on R d and V is spatially 
localized. This scaling was motivated by work of Ilan & Weinstein |16) on the bifurcation of 
nonlinear bound states from continuous spectra for the nonlinear Schrodinger / Gross-Pitaevskii 
equation. The works [161 115] employ the Lyapunov- Schmidt reduction procedure of the present 
work; see also [23117115]. 

The above results describe perturbations of spectra in the case where the perturbing potential 
tends to zero in the strong sense. What of the case where the perturbing potential converges 
weakly? This corresponds to the question of the effective behavior of high-contrast microstructures. 
In [8], the authors consider a class of problems, depending on a small parameter, e, including the 
case where the potential, q^ e '(x) = q(x, x/e), converges weakly as e tends to zero. In particular, we 
considered the small e limit of the spectral, scattering and time-evolution properties for operators 
of the form H^' = — d 2 + q(x, x/e), for a large class of potentials for which y h-> q(-,y) is oscillatory 
(including periodic and certain almost periodic cases) and x h-> q(x, •) is spatially localized. An 
important subtlety arises in the case where q a v(x) = JL q(x, y)dy = 0, i.e. q £ tends to zero weakly. 
In this case, classical homogenization theory breaks down at low energies. Indeed, the homogenized 
operator, obtained by averaging the potential over its fast variations, is Hq = —d 2 , which does 
not capture key spectral and scattering information. Among these are the low energy behavior 
of the transmission coefficient (related to the spectral measure) and the existence of a bifurcating 
bound state at a very small negative energy. We show that the correct behavior is captured by 
an effective Hamiltonian with effective potential well: = —dy -— e 2 A c g(y), A c g(y) > 0. 



Using Theorem 3.1 and the results of [27] . we conclude that has a bound state with negative 
energy of the order e 4 , with a precise expansion for e small. This analysis is then used to derive 
a multi-scale local energy time-decay estimate, for localized initial conditions orthogonal to the 
bound state, in which the different dispersive time-dynamics on different time-scales is explicit. 
These results extend to q^ix) belonging to a class of non-generic potentials (including q^v = 0), 
characterized by the property that q^(x) has a zero energy resonance; see [8] and citations therein. 
The case where q<w(x) is a generic (therefore non-zero) localized potential is treated in 



1.2 Outline, remarks on the proof and future directions 

In section [2] we present background material concerning spectral properties of Schrodinger opera- 
tors with periodic potentials defined on R. In section[3]we give precise technical statements of our 
main results: Theorem 13. II and Theorem 13.41 

The strategy of our proof is to transform the eigenvalue problem, using the appropriate spectral 
transform (Fourier or Floquet-Bloch), to a formulation in frequency (momentum) space. Antici- 
pating a bifurcation from the spectral edge, we break the problem into coupled equations for the 
frequency components located near the band edge and those which are far from the band edge. The 
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precise cutoff depends on the small parameter, A. We employ a Lyapunov-Schmidt type approach 
in which we solve for the far components as a functional of the near components and thereby 
obtain a (nonlocal) reduction to a bifurcation equation for the near- frequency components. Unlike 
classical Lyapunov-Schmidt applications |22| . our reduced equation is infinite dimensional. For A 
small, in an appropriate scaled limit, the bifurcation equation is asymptotically exactly solvable. 
In section [4j we prove a general technical lemma, crucial to the analyses of sections [5] and [6j 
covering the kinds of bifurcation equations which arise. 

Finally, there are several appendices. In appendices [A] and |Bj we present results on regularity 
and other properties of the Floquet-Bloch eigensolutions E^{k), Ub{x\k) = e 2 ™ kx pi,(x; k); see 
also |25| . In appendixjc] we give proofs, by a bootstrap method, of Corollary 3.3 and Corollary 3.6 



which contain more detailed expansions and sharper error terms for the bifurcating eigenstates 
than those in Theorem 13.11 and Theorem 13.41 

We conclude this section with several possible extensions of the present work. 

1. We believe that our methods can be extended to give detailed properties of the spectral 
measure of —d 2 + Q + XV near the band edges. Such information could be used to derive 
the detailed dispersive time-decay, of the sort proved in |5j, but in the present case, for 
initial conditions, which are spectrally localized near band edges. Initial data with spectral 
components away from the band edge can sample a regime where the dispersion relation has 
higher degeneracy, yielding different (slower) dispersive time-decay [3]. 

2. An interesting direction is the extension of the above-discussed results for potentials, q £ , 
which converge weakly to a non-generic localized potential, e.g. g av = 0, to the case where 
q £ converges weakly to a nontrivial periodic potential, Q(x). 

3. Finally, another interesting direction is the extension of the methods of the current paper to 
the study of bifurcations of eigenvalues for multiplicatively small (see |27| ) or weakly conver- 
gent spatially localized perturbations of the two-dimensional periodic Schrodinger operator. 



1.3 Definitions and notation 

We denote by C a constant, which does not depend on the small parameter, A. It may depend on 
norms of Q(x) and V(x), which are assumed finite. C(Ci>C2, • • • ) is a constant depending on the 

parameters £i, ( 2 , We write A < B if A < C B, and A w B if A < B and B < A. 

The methods of this paper employ spectral localization relative to the background operator 
—d 2 + Q(x), where Q(x) is 1— periodic. For the case, Q = 0, we use the classical Fourier transform 
and for Q(x) a non-trivial periodic potential, we use the spectral decomposition of L 2 (R) in terms 
of Floquet-Bloch states; see section [l] and section [2] below. The notations and conventions we use 
are similar to those used in |15j . 

1. For f,g£ I/ 2 (R), the Fourier transform and its inverse are given by 

HfKO = f(0 = f e- 2 ^f(x)dx, F- X {g}{x) = g(x) = f e 2 ^g{i)d^. 
Jr Jm 

2. x an d X are the characteristic functions defined by 

x 4 (o=x(iei<*) = { J; f^tl , xsio = m\< s) = i - xdfi < s) . (i.i5) 

3. T and T -1 denote the Gelfand-Bloch transform and its inverse; see section [2] 

4. L P ' S (R) is the space of functions F : K -> E such that (1 + | • | 2 ) s/2 F G L P (R), endowed with 
the norm 

||F|| L p, s(M) = ||(l + |-| 2 ) a/2 f|| iP(M) <oo, l<p<oo. (1.16) 
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5. W k ' p (R) is the space of functions F : R -> M such that <9jF G L P (R) for < j < k, endowed 
with the norm 

k 

\\ f \\w*.p(R) = E II^H^M < °°> 1 < P < oo- 
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2 Mathematical preliminaries 

In this section we summarize basic results on the spectral theory of Schrodinger operators with 
periodic potentials defined on K. For a detailed discussion, see for example, [101 1251 I2T] , 

2.1 Floquet-Bloch states 

We seek solutions of the k~ pseudo-periodic eigenvalue problem 

(-dl + Q(x))u(x) = Eu(x) , u(x + 1) = e 2mk u(x) , (2.1) 

for k G (—1/2,1/2], the Brillouin zone. Setting u(x;k) = e 2n ' kx p(x; k), we equivalently seek 
eigenfunction solutions of the periodic elliptic boundary value problem: 

(-(d x + 2-Kik) 2 + Q{x))p{x;k) = E{k)p{x;k), p(x + 1; k) = p(x; k) (2.2) 

for each k G (-1/2,1/2]. 

The eigenvalue problem (2.2 1 has a discrete set of eigenpairs {pt(x; k), E b (k)}b>o satisfying ( 1.5 1 
and ( |1.6| |. The functions pb(x;k) may be chosen so that {pb(x; k)}b>o is, for each fixed k G 
(—1/2, 1/2], a complete orthonormal set in Lp Cr ([0, 1]). It can be shown that the set of Floquet- 
Bloch states {u b (x; k) = e 27Tihx p b (x; k), b G N, -1/2 < k < 1/2} is complete in L 2 (R), i.e. for any 
/ G L 2 (M), 

v-^ f 1/2 

fix) - Z_l / ( u b{-,k),f) L 2 m u b (x;k) dk^> 

in L 2 (R) as N f oo. 

Recall further that the spectrum of —d 2 + Q{x) is continuous, and equal to the union of the 



closed intervals, the spectral bands; see (1.7|, 1 1.8 1 



Definition 2.1. We say there is a spectral gap between the b th and (b + l) st bands if 

sup \E b (k)\ < inf |iS&+i(ft)| . 

|fc|<l/2 |fc|<l/2 

Our study of eigenvalue bifurcation from the band edge _E* = Eb„(k*) into a spectral gap, 
requires regularity and detailed properties of Eb(k) near its edges. These are summarized in the 
following results (see a sketch of Eb(k) in Figure [lj left panel). Proofs and references to proofs are 
given in Appendices |A| and [B| for a detailed discussion see [101 1211 ITTI I25| . 

Lemma 2.2. Assume Eb(k t ) is an endpoint of a spectral band of —d 2 + Q(x), which borders on 
a spectral gap. Then fe» G {0, 1/2} and the following results hold: 



G 



1. b even: Eb(0) corresponds to the left (lowermost) end point of the band, 

E b (l/2) corresponds to the right (uppermost) end point, 
b odd: Eb(0) corresponds to the right (uppermost) end point of the band, 
E b (l/2) corresponds to the left (lowermost) end point. 

2. d k E b (k*) = 0; 

3. b even: d%E b (0) > 0, d%E b (l/2) < 0; 
b odd: dlE b [0) < 0, d%E b (l/2) > 0; 

4. dlE b (k,) = 0; 



5. E b (k,) is a simple eigenvalue of the eigenvalue problem (2.11. 

Lemma 2.3. Let E b (ki) denote a simple eigenvalue; thus fci = fc* as above applies. Then, there 
are analytic mappings k h-» E b (k), k H > u b (x; k), with u b normalized, defined for k in a sufficiently 
small complex neighborhood of k\. Moreover, for k real and in this neighborhood (Eb(k),Ub(x; k)) 
are Floquet-Bloch eigenpairs. 

2.2 The Gelfand-Bloch transform 

Given / 6 L 2 (R), we introduce the transform T and its inverse as follows 

T{/0)} = f(x; k) = V e 2 ™*f(k + n), T _1 {/(a;; >)}(x) = e 2 ™ k f(x; k)dk. 

J -1/2 



One can check that 7 1 T = Id. Let 

u(x; k) — e 27Vlkx p(x; k) (2-3) 
denote a Floquet-Bloch mode. Then, by the Poisson summation formula, we have that 
(u(;k),f) L 2 (R) = (p(-,fc),7(-;fc)) i 2 cr([0>1]) . 

Define 

Tb{f}(k) = { Pb (-,k)J(-,k)} L 2 ( [o,i])=/ PbWk)f(x;k)dx. (2.4) 

Jo 

We recall that by completeness of the {p b (x; k)} b >o, the spectral decomposition of / g L 2 (R) in 
terms of Floquet-Bloch states is 

/(*;*) = ^T 6 {/}(%(s;fc) , f(x) = T b {f}(k)u b (x;k)dk. 

b>0 b>0 X / 2 

Recall the Sobolev space, H s , the space of functions with square integrable derivatives up to 
order s. It is natural to construct the following X s norm in terms of Floquet-Bloch states: 



I' ^ JZ^+lbfY \%{<t>}(k)\ 2 dk. (2.5) 

7 -!/2 b>0 



Proposition 2.4. H S (]S.) is isomorphic to X s for s > 0. Moreover, there exist positive constants 
C\, Ca suc/i f/iaf /or aZZ G ff s (R), ?xie fta«e Ci||0||h s (r) < ||0||a- s < ||^||h«{e)- 
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Proof. Since -Bo(O) = inf spec(— d 2 + Q), then Lq = —d 2 + Q — -Eo(O) is a non-negative operator 
and H S (M.) has the equivalent norm defined by ||<^||ir s ~ + Lo) 3 ' 2 (f>\\ L 2. Using orthogonality, it 
follows that 

uw 2 H s ^\\(i+L r /2 n 2 L 2 = J2 f 1/2 \r b {<p}(k)\ 2 \i +E b (k)- E (o)\ s dk 



1/2 



^->n J -1/2 



l{(k}(k)\ 2 dk = \\4>\\%. 



The approximation in the last line follows from the Weyl asymptotics Eb(k) ~ b 2 ; see, for exam- 
ple, [3]. This completes the proof of Proposition 



2.4 



□ 



We conclude this section with a Lemma, which gives various estimates on the Floquet-Bloch 
states of Hq and the spectrum of Hq + w, for a class of periodic potentials, Q, and localized 
potentials, V. These estimates are used within the proof of Theorem |3.4| in section [6] 

Lemma 2.5. Assume that Q G L°° is 1-periodic, and V is such that (1 + | • |)V(-) G L . Let £1 be 
a small neighborhood of hi a simple eigenvalue, such that Lemma\2.3\ applies. Then one has: 



(a) 



feet- 



sup llpbO; fc)|| < sup Y] 



£ 2 ([0,1]) 



< CO, 



(2.6a) 



sup \ \dkPb{-;k)\ 

feeo 



< sup 

feen 



d k Pb(-;k),e 



£ 2 ([0,1]) 



< 00. 



(2.6b) 



sup 



\pb'(x\k)\ \xV(x)\ dx + sup 



\p b ,{x;k)\\V{x)\ dx <oo. (2.6c) 



Proof. We begin by proving that pb(x;k) is uniformly bounded for x G R and k G (—1/2,1/2]. 
Since pb(-;k) is 1— periodic, it will be bounded if its Fourier coefficients are summable. Thus we 
study 

r I 

-2-jvinx 



dx 



^2\{Pb(-,k),e 27Ttn ') L 2 (l0A]) \ = E / P b ( x ' k ) e 

n£Z n€Z J 

Since k G ( — 1/2, 1/2], we can use integration by parts for the Cauchy-Schwarz inequality 

and eqn (|2.2[) for Pb(x; k) to obtain 



L 2 ([0,l])inii 2 ([0,l]) 



\q(x) - Eb(k))p b (x; k) f * — V , </,• 



2^ | <Pb(-; fc),e 2 ™') L 2 ([oa]) | < ||p 6 (a;;fc)| 

ngZ 

+ E 

nez\{o} 



<!+ E i^ll«(0-^(*))p»(-;*)| 



n£Z\{0} 



L2([0,l])' 



Thus, sup fe€( _ 1/2>1/2] ||p6(-;fe)|| i0 o < s up fee( -i/2,i/2]E„ G z |<P6(-;fe),e 2,r " l ')L2([o,i])| < oo. 

Let us now turn to the study of d k Pb(x; k) in (b). Deriving \2.2\ with respect to the parameter 
k yields 



(d x + 2nik) 2 + Q{x)) d kPb (x; k) = E b {k)d k p b {x- k) + (d k E b (k) + 4md x - 8n 2 k)p b (x; k). 



<S 



Following the same method as above, one obtains 

||<9fcP(,(-;fc)|| LOO < ^2 \i d kPb(-; k), e 2 "™') L 2 ([0jl]) | 

< C(\\Q\\ L ^,E b (k))\\d kPb (x; k)\\ L2([01]) + C(d k E b (k)). 

The finiteness of \\dk P b(x; k)\\ 2 and dkEt(k) for k G f2 is a consequence of Lemma 
(b) follows. Finally, (c) is a straightforward consequence of (a), with (1 + | • G L 1 . 



2.3 



thus 

□ 



3 Bifurcation of defect states into gaps; main results 

Consider the eigenvalue problem: (-<9 2 + Q(x) + XV(x)) ip x = E x ip x , V S L 2 (R), where Q{x) 
is 1— periodic, A is small, and V(x) is spatially localized. Our first result concerns that case where 
Q = 0: 

Theorem 3.1 (Q = 0). Let V be such that V G W 1 '°°(R); thus f R (l + \x\)\V (x)\ dx < oo suffices. 
Assume V(0) — f R V < 0. There exist positive constants Ao and C(V, Ao), such that for all 
< A < Ao, there exists an eigenpair (E x ,ip x ), solution of the eigenvalue problem 

{-dt + \V{x))ip x {x) = E x iP x (x) (3.1) 

with negative eigenvalue of the order A 2 . Specifically, 

\ 2 

[ ~ 



E 



4 



' CA 5/2 



sup 



ip x (x) 



-PI \ 



V\\x\\ ' CA 1/2 . 



(3.2) 
(3.3) 



The eigenvalue, E , is unique in the neighborhood defined by (3.2|, and the corresponding eigen 
function, ip, is unique up to a multiplicative constant. 



Remark 3.2. Theorem \3.1\ shows, and is essentially proved by demonstrating, that for small 
positive A, the leading order behavior of the eigenstate (E x ,ip x (x)) is a scaling of the unique 
eigenstate of the attractive Dirac delta potential: 

( E x , V A M ) w ( A 2 #o , So(Az) ) , 

where 9 = V > and g (y) = e~ e ° M satisfy 



-K + 



V ■ S(y) 



9o{y) = -d 9o(y)- 



(3.4) 



The error bounds in Theorem |3.1| are not optimal. However, the bootstrap argument of Ap- 
pendix [c] can be used to recover a higher order expansion on E x , similar to that obtained in |27j . 

Corollary 3.3. Assume (1 + \x\ 2 )V G L 1 , and V{0) = J* B V(z) dz < . Then E x , as defined 
in Theorem\3.1\ satisfies the precise estimate: 



E 



-A 2 [61(A)] 2 , with9(X) = -\ I V -\\ f[ V{x)\x-y\V(y) dxdy + £>(A 3/2 ) . (3.5) 

2 J M 4 JJ R 2 
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Simon |27] and Klaus |18| prove expansion (3.51, under the conditions: (1 + |a;|)|V / (a;) G L (R) 
and L V < 0, with the error term o(A). Corollary |3.3| gives a sharper error term under a more 
stringent decay condition on V. The proof of Corollary |3,3| is given in Appendix [C] . 



Theorem 3.4 ( Q non-trivial, 1-periodic). 
(1 + | ■ |)V(-) G L 1 and V G IT. Let E K 



Let Q G L°° be 1-periodic, and let V be such that 
: k G (—1/2,1/2] — > R denote the band dispersion 



function associated with the band of the continuous spectrum of — + Q(x). Fix a sp & 
band edge of the (6*) th band; thus E+ = E bt (k*), where fc* = or fc* = 1/2 (see Lemma 2.2). 
Assume either 



d 2 k E bt (k t )>Q and 



r 

J — ( 



\ub,{x;k„)\ 2 V{x)dx < 0, 



(3.6) 



d 2 k E bt {k t )<Q and 



\u K (x\k,)\ 2 V(x)dx > 0. 



(3.7) 



Then, there exist positive constants, Ao and C — C(Ao, V, Q), such that for all A < Ao, the following 
assertions hold: 



1. There exists an eigenpair ( y E x ,tp (cc)) of the eigenvalue problem 

(-d 2 +Q(x) + XV(x))^ x (x) = E x i^ x (x), tp x &L 2 (R) 

2. E x and ip x (x) satisfy the following approximations: 

\e x - ( E K (k,) + X 2 E 2 ) j < CA 2+1/4 . 
sup 4> X (x) — Ub, (x; fc*) exp(AaoH) < CA 1,/4 . 



where 



E-, 



\jr oo \u bt (x;k t )\ 2 V(x)da 
2^dlE b ,(k t ) 



fZ \u b .(x;k.)\ 2 V(x)da 
^dlE b „(k*) 



(3.8) 

(3.9) 
(3.10) 

(3.11) 



(3.61 and (3.71, we have that ao < 0. Note also that the direction of bifurcation of E is 
sgn(£ 2 )= -sgn(dlE b ,(k*)) . 



given by 



3. The eigenvalue, E , is unique in the neighborhood defined by (3.91, and the corresponding 
eigenf unction, ip x , is unique up to a multiplicative constant. 



Remark 3.5. By Theorem 3.4 the bifurcating eigenvalue E x lies in the spectral gap of —d 2 + Q(x) 
at a distance 0{\ 2 ) near the spectral edge E„; see Figure^ Moreover, E2 is the unique eigenvalue 
and g {y) = e a ° M is the unique (up to multiplication by a constant) eigenfunction of the effective 
(homogenized) Hamiltonian: 



d 

dy 87r 



\d 2 k E K {k*) 4- + 



dy 



r 



\ub m (x; k*)\ V(x)dxx5(y). 



The following refinement of Theorem |3 . 4| can be proved via the bootstrap argument presented 
in Appendix [C] 
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Corollary 3.6. Assume (1 + |a;| 2 )V £ L 1 , and hypotheses of Theorem ( 3.4 1 hold. Then one has 



E x -E b ,(k,) = \ 2 (E 2 + \E 3 ) + 0(X 3+1/4 )= -A 2 a ^T (M [9(A)] 2 , (3.12) 
where Eo is as in i |3~TTj ), 

e 3 = ( ^ gi; 8 ^ ))2 (/°° K fc*)! 2 v(x) dx) v(aOK to Ml> - 2/IK (2/; 



and 



1 



9(A) =-- K* {x; K)\'V{x) dx (3.13) 

-\ X ffi^fi. \ II V(x)\u bt (x;k t )\ 2 \x - y\\u bt (y;k,)\ 2 V(y) dx dy + 0(X 1+1/4 ). 
4 d£E b ,{k*) JJ R 2 

Remark 3.7. For the case Q = 0, the spectrum consists of only one semi-infinite band which we 
can label the b — band. In this case, uo(x; A;* = 0) = 1 and Eo(k) = 4-7r 2 fc 2 . Therefore, to leading 



order, relation (3.131 simplifies to the result of Corollary 3.3 and the two results are consistent 



4 Key general technical results 

In this section, we study the operator Co[9], defined by: 

7(0 h+ C [d]f(0 = (4^ 2 ^ 2 + 6 2 ) /(O -By (\t\ < \- p )J^x(\v\ < A^) f(v) dq. (4.1) 

Here, A, B and (3 are fixed positive constants. The operator £-o[0] appears in the bifurcation 
equations we derived via the Lyapunov-Schmidt reduction; see section |l~2| 
In space, we have that Co[6] is a rank one perturbation of —Ad y + 2 : 

CoW EE (-Ad 2 + e 2 )f(y) - g <^ sinc f^L .) sinc , (4 .2) 

where sinc(;z) = sin(2)/z. Cq[8] is a regularization of the operator: 

( h a ' b + e 2 ) f = {~Ad 2 y - BS( y ) + e 2 ) f , (4.3) 

appearing in the effective equations governing the leading order behavior of bifurcating eigenstates; 
see Remarks 13.21 and 13.51 



4.1 The operator £ 

Lemma 4.1. Fix constants A > 0, B > and (3 > 0. Define, for 9 2 > 0, i/ie linear operator 
f{t) ^ C [8]f(0 = (4n 2 Af + 6 2 ) f(0 - B X (jf | < A"' 3 ) J X (\v\ < A"") /(*?) *j. (4.4) 



iVote t/iai £ [6>] : l/QR) -> L 1 '" 2 (R); see Og| 



_Z. There exists a unique 9 > smc/i i/iat £o[#o] ftos a non-trivial kernel. 
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2. The "eigenvalue" 9 is the unique positive solution of 

x(|g|<A- g ) 

4tT 2 A£ 2 + 9' 2 



B 



3. The kernel of £o[8o] is given by: 

kernel = span 



{f (0}, where / (£) 



x(kl<A-Q 



4- 6o = #o(A) can be approximated as follows: 



do (A) 

5. Define g(x) = exp(ao|a;|), with cto = 

sup J r_1 {/o} (a;) 



B 



2^/A 



o-T < 0. Then one has 



-g(x) | <C(A ) B)A* 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



Proof. First note, by rearranging terms in the equation £o[#o]<? = 0, that any element, <?(£), of the 
kernel of £ [6»], is a constant multiple of the function / A (£; 9) = < A _,s ) x (4tt 2 y4£ 2 + 6> 2 )~\ 
Thus, if g is non-trivial then it is strictly positive or strictly negative and therefore f R g ^ 0. Next, 
note that a necessary condition for g to lie in the kernel of Co[9] is that equation ( |4.5[ ) holds. To 
see this, divide the equation £o[#o]ff = by (47r 2 ^4£ 2 + 9q 2 ) and integrate d£ over K. This yields: 



r° * x [ i - b r 

J — OO •/ — ( 



x(iei<A^) 
4-ir 2 A£ 2 + 6 2 



= 0, 



(4.9) 



By the above discussion, if ~g is non-trivial then J R g ^ 0. Hence 9 2 satisfies J{9 2 
B fZo = 0. Since J : (0,oo) - 



R is smooth, J'(X) > 0, lim J(X) = — oo and 
x->o 

lim J(X) = 1, the function J has a unique positive root, which we denote by (J 2 ,. One can check 
by direct substitution and the condition J(9 2 ) — 0, that any multiple of 

/o(0 = A(£;0 O ) = x(l£l < A"") x {A^Ae+91)- 1 

satisfies jCo[#o]/o(£) ~ . 



(4.10) 



The approximation to #o(A), (4.7|, is obtained as follows. Let 9q denote the unique solution of 
J(0q) — and 9g its positive square root. Then, 



1 = / 

B Jm 



4tt 2 i^ 2 +^ ? 



/ 

J3t 



l + (x(KI<A-")-l) 



4tt 2 A £ 2 + 6> 2 



2\/]4 6*o 



The last term can be bounded as follows: 

i-x(iei<A-*) 



4tt 2 A£ 2 + 91 



ci>a-p 



4tt 2 A £ 2 + (9 2 



+ 



X(|^| <A^)-1 
47rM£ 2 + 

(4.11) 



t\>\-f> 



A* 3 



4tt 2 A £ 2 ~ 2n 2 A 



(4.12) 



Relations (4.111, (4.121, after rearrangement of terms, yield ( 4. T I . 
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Finally, let us turn to the asymptotic expression for T 1 |/o| {%) given in (4.8 1. By residue 

— ao _ B/2 



computation, one has 

2 



47r 2 |<e| 2 + a§ 



It follows that 



sup T 



< 



1 {fo} (x) 



B 



9(x) < /o 



2 _ 
B 9 



x(\t\<x- p ) 



X(ICI < A- 



^ 2 A\tf + f^ 



4ttM£ 2 + 6» 2 



^ 2 m 2 + S 



4n 2 A£ 2 + 9 2 

l-x(ICI<A- fi 



4tt 2 ^ 2 + g 



A bound on the second term follows from (4.12 I. The first term is easily bounded, using (4.7 1, by 
C(A,B)\ fi 1 for some constant C(A,B) > 0. Estimate ( |4.8| l follows, and the proof of Lemma 
is now complete. 



4.1 



□ 

We shall also require a result on the solvability of the inhomogeneous equation 

(£o[0o]£) (£)=%), (4.13) 



where £o[#o] is defined in (4.4|. 

Lemma 4.2. Let h be such that x (ICI < A _/3 ) = and satisfy the orthogonality condition 



fo,h) = 0, 



(4.14) 



where fo, displayed in (4.61, spans the kernel of C[9q]. Then, 



1. any solution of the inhomogeneous equation (4.131 is of the form 

m = (c+h(o)Mo = (c+h(o) x + [I , (4.i5) 

for some constant C 

2. The unique solution of ( |4.13[ ) such that J R = is obtained by choosing C — 0: 

m = %(o M)- (4.i6) 

Proof. To show that ( |4.15| l solves the inhomogeneous equation ( |4.13[ ) we simply insert the func- 



tion (4.151 into (4.131, and use the properties: C{9o)fo = and ^/o, hj 2 = 0. This gives 

(£o[9o]$) (0 = {^ 2 A( 2 + 6' 2 )h(Ofo(0 - B x (ill < A"- 9 ) J°° h(v)fo(v)dr] 

x(KI<A-")h(0 



(4tt^ + 0$) 



4tt 2 A£ 2 + 9 2 



b x iei<A 



fa,h 



HO- 



The converse clearly holds by Lemma |4.1| since the difference of solutions of the inhomoge- 
neous equation solves the homogeneous equation (4.4 1. Finally, using the orthogonality condition 

\fo,h) =0, one has that J R cp = C J R fo = if and only if C = 0. □ 
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4.2 A perturbation result for £ 

As discussed in the introduction, our strategy is to obtain a reduction of the eigenvalue problem 
for Hq+w to an eigenvalue problem (the bifurcation equation) for functions supported at energies 
near the band-edge. These reduced equations have a general form which we study in this section. 
Let Zi and Z2 denote Banach spaces with Zi, Z% C L\ oc . Assume that for any (/, g) G Z\ x £2, 



\(f,g)i 



\f9\\ 



ll/l 



z 2 



and 



(i + cTVIU^H/IU- 

(4.17) 

Furthermore, we also assume that /o G Z\ n Z2, where ( 8$, fo] is the unique normalized solution 

of the homogeneous equation Ca[0]f = 0; see Lemma 
We seek a solution of the equation: 



4.1 



£o{e\f = R{f) 



(4.18) 



where Cq{9) is the operator defined in ( 4.4 1 and the mapping / i-> R(f) is linear and satisfies the 
following properties: 
Assumptions R a ,/3: 

There exist constants a > 0, ft > and Cr > such that for any / G Z2 



X |£|<A^ )R{f)(0 = R(f)(0 



and 



•Zl 



< C R X a \\f 



\z 2 



(4.19) 



where 



In the above setting we have the following 

Lemma 4.3. Let (0o, fo{£,)) be the solution of £o(#o)/o = 0, as defined in Lemma 4-1 
A, B and ft > are fixed. Let R : f £ Z2 — ¥ Z\ be a linear mapping satisfying Assumptions 
Ra,p displayed in (4.191, where Z\,Z2 satisfy (4.171. Then there exists Xq > such that for any 
< A < Ao, the following holds: 



1. There exists a unique solution y0,f(£) J G R x Z2 of the equation (4.181, such that 

/oo 
no -ho at 
-00 



0. 



< CA° 



with C = C(A, B,Cr, ft), independent of X. 
2. Moreover, one has /(£) = \ < T(0 an d \0 2 

Remark 4.4. To prove Theorems \3.1\ and \3.4\ we shall apply Lemma \4-3\ to the operators: 

1. -dl + XV{x) with (Z!,Z 2 ) = (I/ 00 , L 1 ), and 

2. —dx + Q(x) + XV(x) with (Z\,Z2) = (L 2 ' _1 , L 2 ' 1 ), where L 2 ' s is the space of locally integrable 

||( 1 + |£|2y./2l 



functions such that 



\F 



L2- 



It is straightforward to check that such spaces satisfy ( |4.17[ ), and fo G Z\ n Z2- 

Proof of Lemma \4-3\ Our strategy is to use a fixed point argument. We seek a solution (9 2 ,f ) 
to (14. 181) of the form 



and 



f = fo + fi 
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Any solution / of filty satisfies /(£ ) = % < A _/3 ) /(£)• Therefore, since / (f ) = X < A _/3 ) /o 
it follows that /i(£) = X (ICI < A* 3 ) Substitution of these expressions into ( 4.18[ ) yields 



(4^^ + ^) X (|f|<A-' 3 ) (/o + /i)«) 



X |€| < A" 



B 



X(N <A" /3 ) (/o + /i)(r;)dr ? = J R(/o + /i)(0 



Rearranging terms yields the following equation for /i , in which Q\ is a parameter to be determined: 
(£>M/i) (0 = -0? (/o + /i) (0 +R(fo + fx) (fl. (4-20) 



By Lemma 4.2 (4.201 is solvable in L only if the right hand side is L - orthogonal to /o: 



h,-el [fo + fi)+R(f + h 



L 2 



= 0. 



Solving for 9i , we obtain 



L 2 



(4.21) 



In summary, equation (4.181 can be rewritten equivalently as two coupled equations in terms of 

fi and 6f: |4!20| -( [42T| . 

Substitution of 9 1 in (4.201, or equivalently projecting the right hand side of (4.201 onto the 
orthogonal complement of span{/o}, yields the following closed equation for /i: 

' To, R( fo + fx) 



(£ o [0o]/i) (0 



L- 



foJo) L2 + {fo,h /L 



(fo + fi)(0 + R(fo + fi)(0- (4-22) 



By Lemma 4.2 /i is a solution of ( |4.22[ ) with J R fi — if and only if : 

A(0 = £f(A)(0, 

where 



(4.23) 



X(I^I<A-Q 

47T 2 Af 2 + 6» ( 2 



V 



fo,R (fo + fx)) 



fo + fx )(t) + R (fo + fx (0 • (4.24) 



We solve the fixed point equation (|4.23[ ) by the contraction mapping principle. Once /i has been 
obtained, (?! is determined using ( |4.21[ ). 
Introduce 



e z 2 



11/1 



32 



< C«A 



for some fixed Ch > 



(4-25) 



Note that S is a closed subset of the Banach space Z-i. We next show that there exists Ao > 
such that for all < A < Ao: G '■ S — > S and Q is a contraction mapping. As a consequence, it 
will follow that for < A < Ao, there is a unique solution /i £ S of the equation /i = G{fi) and 
therefore of (4.221. Moreover, ||/i|| < A" by definition of S, and one can check: 

/ -/fo,R( fo + fx] 



[fx=[ S{fx) = ( 

JR JR JR 



MO 



L 2 



fo,f ) L2 +{fo,fx /L2 



fo + fx) (0 + R( fo + fx) (0 d£ 
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It then remains to obtain an estimate of 9i 2 = 9q 2 — 9 2 . From (4.21 1, one has 

\el\ < \(fo,R(fo + fi)) I - — 1 ,. < A° 



where we used (4.171 and (4.191, and the fact that for A sufficiently small, <^/o,/o J o > c > 0, 
where c is independent of A. Lemma |4.3| is proved. □ 

Proof that Q : S — > S is a contraction mapping: The result is a consequence of the two following 
claims, proved below: 

Claim 4.5. There exists C H = C (e ,A,C R , ||7o|| z J > such that \\G(0)\\ Z2 < \C H \ a . 



Claim 4.6. There exists A > such that ifO < A < A , then \\Q(fi) - S(/a) < § /i - h\ 



^2 



Indeed, Q maps 5 = -|/ (E Z<z : ||/|| z < ChA q j into S 



< ||e(/)-e(0)|| Z2 + ||g(0)|| Z2 < -||/-0|| Z2 + -ChA" < C„A° 
and Q : S — s> S is a contraction mapping by Claim |4~6| 
Proof of Claim \Jp% By definition, one has 

x(fcl<A-' 



/o,i«(/o)) i2 M)+i«(/o) (0) 



4ttM£ 2 + 6» 2 

It follows, from our assumptions ( |4.17[ ) on functional spaces (Z\ , 22 



ie(o)|| 



22 



X(I£I<A- 8 ) 
4tt 2 ^ 2 + 9\ 



h,R[h)) L2 \\f4 z ,+ 



x \\i\<\-"){i+e) 






4tt 2 A£ 2 + 61 




1 + 1 ■ I 2 
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Claim 4.5 is now obvious, using the smallness hypothesis on the operator R, (4.191 



(4.26) 

□ 



Proof of Claim \4^6\ Let us decompose the mapping Q as follows 
X (ICI < A-' 9 ) / (To,R(fo + fi 



9{h - h 



4tt 2 A£ 2 + 6 2 



fo + fi (0 



fo,R[fo + h 



L- 



(fo + ?*) K) + R (h - h) (£) 



4tt 2 A£ 2 + 6»g 



4tt 2 24£ 2 + 2 
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The following estimate follows from our Assumptions ( 4.17[ ) on the spaces (Z\ , Z; 



0(/i-/a) 



I2 2 



47r 2 A£ 2 + 2 



2v 



47r 2 A£ 2 + 6» 2 



^2 



x(iei<A-^) 



47t 2 a^ 2 + el 



\Si[fi] - S 1 [f 2 ] 



< \\S 1 [f 1 ]-S 1 [f 2 ] \ + \\R [h-h 
II \\z 2 II 



\z 2 



I Si 



x(|£l<A^)(l + £ 2 ) 




R (I - h) 


4tt 2 A£ 2 + 6» 2 




1 + | ■ | 2 



z 2 



The second term in (4.271 is estimated using assumptions R a .p, (4.191: 

||-R (/i - £) |l _ < c R x a \\f l -U ■ 



\z x 



\Z 2 



(4.27) 



(4.28) 



Let us now turn to the first term in (j-4.27 1 
Si[h]-Sx[f 2 } = 



fo +fi) + 



h,R[h + h 



fo,fo) L2+ {fo,h /L ^ 



fo + h 



{fo,R(n-h)) L2 (f +h) {f ,R(f +h)) L2 (n-h) 

( 



L 2 



fo,fo) +(foJ: 



L- 



L- 



fo,R(fo + h)) L2 (/0 + /2) 



1 



/0,/0) +(/0,/2 



= 1 + 11 + III. 

The result is a consequence of the following estimates: 

{fo,g) L2 < C\\f \\ Zi \\g\\ Z2 < dHffll^, 
[fo,R(g)) L2 < C\\f \\ Z2 \\R(g)\\ Zi < C 2 \ a \\g\\ Z2 , 

with Ci = C*(||/o|| Zi ) and C 2 = C2(||/o||^ 2 ,Cb). 

Using the above estimates, one easily checks that, provided CiA Q < 1/2, 

||/|| Za < 2C7 2 A Q ||/ 1 -/ 2 || (\\fo\\ z +C H \ a ), 



(4.29) 



3 2 



\\II\\ Z2 < 2C7 2 A Q (||/ || z +C H \ a ) \\fi-h\ 



z 2 



\\HI\\„ < 4Ci /1-/2 C 2 A Q (||/o|L +C H \° 
2 II }\z 2 "^2 



Thus if CiA Q < 1/2, one has 



Si[/i]-Si[/ 2 ]||_ < IKIU + II^IU + II^IU, < ^\\h-f: 



\z 2 



z 2 



*z 2 



\z 2 



(4.30) 



Plugging ( |4.28[ ) and ( |4.30[ ) into ( |4.27[ ), it follows the existence of a constant, Co > 0, such that 

||<5(/i) — G(f 2 )\\z 2 < CoA a ||/i — f 2 \\z 2 - Thus for < A < Ao < C °, we obtain a contraction and 
claim [476] is proved. □ 
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5 Proof of Th'm 3.1; Edge bifurcations for —d\ + XV (x) 

In this section we study the bifurcation of solutions to the eigenvalue problem 

{-dl + XV (x)) 4> x (x) = E x ip x (x), ^ A 6l 2 (R) (5.1) 

into the interval (— oo,0), the semi-infinite spectral gap of Hq = —d 2 , for V a spatially localized 
potential, and A > sufficiently small. 

More precisely, we prove Theorem |3.1| which is a particular case of our main result, Theo- 
rem 3.4 In this case (Q = 0), the Floquet-Bloch eigenfunctions are explicit exponentials, thus 



calculations are more straightforward and error bounds on the approximations are sharper. 

5.1 Near and far frequency decomposition 

Taking the Fourier transform of ( |5.1[ ) yields 

4n 2 e - E x ) + X [ V (£ - C)^(C) dC, = 0. 



(5.2) 



We shall study (5.2 I via an equivalent system for the 

near frequency component: {tp x (^) : |£| < A r } and 
far frequency component: {?/> A (£) : |£| > A r } of i[) X . 



Let r be a parameter, at first chosen to satisfy: < r < 1. Below, we shall further constrain 



the range of r— values. Recall the cutoff functions, % and Xi introduced in ( 1.15 I and 

xyr (0 = x(x~ r 0, With 1 = Y (0 + v K). 



Multiplying (5.2 1 by this identity we get 



= (4tt 2 |£| 2 - E x ) (x xr (0 + x xr (0) 

+ A /°° (x,M) + x x A0)v(£-0 (xxA0 + x xr (0)^(0dC 

Introduce notation for the near- and far- frequency components of t/j X - 

&«,(0 = x xr (0^(0 and ^r(0 = x x A0^ x (0- (5.3) 

Then, the eigenvalue equation is equivalent to the following coupled system: 

(5.4) 
(5.5) 

In what follows we shall set E x — —X 2 6 2 , where 9 — 0(X) is expected to be 0(1) as A J, 0. This 
anticipates that the bifurcating eigenvalue, E x , will be real, negative and 0(X 2 ). 



Un 2 \£\ 2 - E X ) Vw r (£) + X Xx r(0 / V(C- C)(^ne ar (C)+^far(C)K = 0, 

J — oo 

(47T 2 |C| 2 -i5 A )^ ar (0 + AX Ar (0 ^ %-C)(fnear(C)+^a r (C))rfC = 0. 
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5.2 Analysis of the far frequency components 



We view ( |5.5[ ) as an equation for iptax, depending on "parameters" (i/'near; A). The following propo- 
sition studies the mapping (V>near; A) i-> Vfar. 

Proposition 5.1. £e£ -tpnear £ i 1 . There exists Ao > 0, suc/i t/iai /or < A < Ao, f/ie following 
holds. Set E x = —\ 2 8 2 , with \9\ < 7rA r_1 , r G (0,1) . There exists a unique solution, "tpf ar , 
of the far frequency equation (5.5 1. The mapping (ipnenr',X) >-> ^/«r [V'near; a1 maps L 1 (R) x R to 
L 1 (R) and satisfies the bound: 



W^/orWzA < CQ\V\\loo) A 1 ""- U near \\ L i. (5.6) 

Proof. We seek to solve | ]5.5| for ip{ al . as a functional of Tp ncs . v . Since |£| > A r , with < r < 1, 
and |0| < TrA* --1 , we have^r 2 £ 2 - E x \ = |4tt 2 £ 2 + A 2 6» 2 | > 37r 2 A 2r , which is bounded away from 
zero for any fixed A > 0. Dividing ( |5.5[ l by 47r 2 £ 2 - E x = 47r 2 £ 2 + A 2 6» 2 , we obtain that ( |5.5[ ) is 
equivalent to the equation: 



where 



I + Tx) Vfer(l) = " TA^near (£) , (5.7) 



(Ta g) K) = J Kx(M) 9(0 d( and /C A (£,0 = V^ 2 + A^ ^ ~ • 
We next show that the integral operator 7a, viewed as an operator from L 1 to L has s mall 



norm, for A small. This implies the invertibility of / + 7a and the assertions of Proposition 
Let j e i 1 . One has 

\\fxg\\ Ll < C\j ,_ 9 ^ 1 , 9n9 <% \\V\\ roo \\g\\ rl < X 1 ^ \\V\\ roa \\g\\ rl . 



5.1 



«|>V 4ir 2 £ 2 + \ 2 6 2 



Thus 7a is bounded from L 1 to L 1 with norm bound: HTaII^^x < A 1 "" \\V\\ L , X . For r G (0,1), 
IITaIIl 1 ^!, 1 — > as A — > 0. Therefore J + 7a is invertible, for A sufficiently small. Moreover, 



^near r 1 5 



which implies the bound (|5.6l). Proposition |5. 1| is proved. □ 



5.3 Analysis of the near frequency components 



Now that we have constructed ^>f ar as a functional of linear and A (Proposition 5.1 1, it is possible 



to treat (5.4 1, for A small, as a closed equation for a low frequency projected eigenstate, V'near(Ci A), 
and corresponding eigenvalue E x . Substitution of i/)f ar = ^far [V'near , A] into (5.4 1 yields: 



(4^|C| 2 - E X ) VW(£) + \ Xx r(0 J V{£ - C)Vw(CR + AJCV(0«(0 = ( 5 - § ) 

where R is defined by 

%) = C) $ta$«r,A](C) d( . (5.9) 

Recall that i/>f ar [ifioaut > A] is in L 1 , and of size O ^A 1_r HV'nearlU 1 ^ by Proposition jsTTj 
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Our next goal is, via appropriate expansion, reorganization and scaling, to re-express (5.81 
as a simple leading order asymptotic equation plus controllable corrections. The terms in (5.81 
are supported in the near (low) frequency regime. Note that for |£| < A r and |£| < A r we have 
II - CI < III + ICI < 2A r . Taylor expansion of 9(£ - C) gives 9(£ - C) = ^(0) + (| - C)9'(v), for 
some n = t?(C,0 such that < 2A r . Using this expansion in the second term of (5.81 yields 



where 1Z 



(4tt 2 |£| 2 - £ A ) &«x(0 + Ax^(O^(O)/ ^ ear (C)dC = Axa-(I)^ pnca r ; a] ($, (5.10) 

near ■ 



7ii + 7e 2 , with 



7fe(£) = -%) = -fvit-Q V^ far [^nea r , A] (C) 7k(0 = - | (| - C) V' (l)i, r (C) <*C 



We now introduce the scaled near-frequency Fourier component, $a, by 



*/w(l;A) = ( | 



A 



Note that 



-$A 



A 



= $A 



(5.11) 



(5.12) 



We also denote E x = — \ 2 2 , and restrict to # = 0(X) satisfying the constraint in the hypotheses 
of Proposition 5.1 Substitution of (5.111 into (5.101, defining £' = A£ and dividing by A yields 



(5.13) 



the following rescaled near-frequency equation: 

(^^'l 2 + 9 2 ) 8 A (0 + Xxr-l (e')V(O) ^ $A(C'X = XAr-I (O^'C** )(0 
where = ft [Vw;A] (A£') = 7*i(0 +*a(0. with 

U(Ae' - C) ^farl^ncar, A] (£) d( , (5-14) 

(A|' - C)^'(??)Vw(C) dC = -A^(|' - C')^'W$a(C'K' ■ (5-15) 



74(1' 



Equation (5.13 1 is in the form of the class of equations to which Lemma 4.3 applies. We shall 
use Lemma ■ 



4.3 



to obtain a non-trivial eigenpair solution ( <&\,9(\) ) of (5.131. Toward verification 
of the hypotheses of Lemma 4.3 we next bound the right hand side of (5.131. 

+ IIV* 



Proposition 5.2. Let V be such that V 



W 1 '" 



hand side of the rescaled near-frequency equation ( 5.13[ ) satisfies the bound 

\\xxr-im'fr)\\ LBB < C(\\9\\ wli00 ) (A^ + A") \\$ x \\ Li 

Proof. We proceed by estimating each term individually. 
Estimation oflZ'i(£'), given by (5.14 1: By Proposition |5.l| one has 

i 1— r || T II 



< oo . Then, the right 

id 

(5.16) 



^far 



[</W,A]|| < C(\\V\\ L ~) A 1 



(5.17) 
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Plugging (5.17l into (5.141, and making use of (5.121, we have 



l^il 



F(A£'-C)Vwtyw,A](C) rfC 

< \\V\\t,~ ||^[^ear,A] IL, < C(\\V\\ L ~) 



It 1 ■ 



Estimation o/7?-2(C )> given by (5.151: We have the bound 



1x^(0^1 



Xxr-i(0 HZ'-OV'(ri)$ x (Odt'\\ L ~ < 2A r ||V'l 



using that $a(C') = Xv-i {C)$x(C), so that l£' - CI < 2A r ~ 1 . Proposition 5.2 is proved. □ 
Remark 5.3. We expect that by using a higher order Taylor approximation of V(£ — in the 



second term of equation (5.8 1, it should be possible to obtain a variant of Proposition 



bound which is higher order in X. This would require a higher order variant of Lemma 

5.4 Completion of the proof 



5.2 with a 



4.3 



We now prove Theorem |3.1| by an application of Lemma 4.3 to equation ( 5.13 1, using the remainder 
estimate of Proposition |5.2| 

Proof of Theorem\3.1\ We construct ip X , solution to (|5.2[) as ijj x — -ip{ &r + ^ ncar , where V>f ar , ^near 



satisfy ( 5.4 H 5.5 1. The far-fr eque ncy component, ipc^, is uniquely determined by ^ near and A suffi- 



ciently small; see Proposition 5.1 Now set V) ne ar(0 = y ^a (j )• Since V £ W^ 1 ' 00 , Proposition 
implies that the rescaled near-frequency equation ( 5.13[ l can be written as 



5.2 



(4^|£'| 2 + 9 2 ) *a(0 + Xxr-i(t')V(0) J <&a(CK = Xxr-t(?)7l(* x )((!), 



(5.18) 



with ||7?-(w)|| £ oo < C X a \\u\\ L i, where a — min(l — r,r) and C = C(|| V|| w i,oo). From now on, 
we set 

r = 1/2 = a 



as this yields optimal estimates. Applying Lemma 4.3 to (5.181 with 



A = 1 and — B = V(0) = L V (assumed to be negative), we deduce that there exists a solution 
^ 2 ,$a^ of the rescaled near- frequency equation ( |5.18[ ), satisfying 



||$A-/o||x,i < CA5 and \9 2 - 9%\ < C A 5 . 
Here ^#o(A),/oj is the unique (normalized) solution of the homogeneous equation 

£o,a (Oo, fo) = (4^ 2 + 9 2 )fo + X (^1 <A-*)V(0) fx(\v\ <A-i)/ (r ? )d7 ? = 0, 



(5.19) 



as described in Lemma 



4.1 



Thus V'near(C) = i5 A (§)and£: A = ~X 2 9 2 (X) are well-defined 

A ,/.A\ 



In conclusion, the eigenpair solution to (|5.2[) (i.e. (|3.1[)) , (£7 ), is uniquely determined by 



-X 9 (X), and V = J 7- 

near + V'far). 



21 



Estimate (3.2 I, the small A expansion of the eigenvalue E , follows from 15.191. The approx- 
imation, (3.3 1, of the corresponding eigenstate, tp x = i/> nC ar + t/'far, is obtained as follows. First, 
by ( |5.19[ ) we have 



■0near(?7) - A 



4TT 2 \rj\ 2 + \ 2 8$ 



M 



< A 1/2 . 



The high frequency components are small, as is seen from the following calculation: 



< A 



dr) 



„|>Al/2 4^21^1- 



47r 2 |r7| 2 + A 2 V(0) 2 
Finally, from Proposition |5.1| one has (with r = 1/2) 

and | linear 1 1 L i = ||$a | Li -» 11/o Li (as A -» 0). 
Altogether, ( |5.20[ ), ( |5.21| ) and ( |5.22[ ) yield 



< A 1/2 . 



J"" 1 A 



47T 2 



+ A 2 2 



V; A - A 



47T 2 | 



A 2 2 



(5.20) 



(5.21) 



(5.22) 



< A 1/2 . 



Note, by residue computation, that T 1 {(47r 2 j • | 2 + A 2 0, 2 ) 1 } = (A6» ) 1 exp(— Aflojd), with 



| L V > . Thus estimate (3.31 holds. This completes the proof of Theorem 



3.1 



□ 



6 Proof of Th'm ^ Edge bifurcations of -d* + Q + XV 

Let Q(x) denote a non-trivial 1-periodic function, Q(x + 1) = Q(x). In this section we study the 
bifurcation of solutions to the eigenvalue problem 



{-d 2 x + Q(x) + XV(x)) tp*(x) = E x ip x (x), V G L 2 (R) 



(6.1) 



into the spectral gaps of — <9 2 + Q(x). We shall proceed by the same general approach of section [5] 
That is, by appropriate spectral localization, in this case by applying the Floquet-Bloch transform, 
we reduce (6.1 1 to an equivalent near-frequency eigenvalue problem supported on frequencies lying 
near a spectral band edge of — <9 2 + Q(x). 



6.1 Near and far frequency components 

We take the Gelfand-Bloch transform of ( |6.1[ ) and get 

- (d x + 2irik) 2 ^(x; k) + Q(x)^ x (x; k) + A (Vip)~ (x; k) = E X ^(x; k), (6.2) 

where 

(VV> A ) ~ (x; k) = J2 ( V ^) A ( k + n ) = J2 e2 "' nX iy *^ X )(k + n). 

Here, the quasi-momentum, k, varies over the interval (—1/2, 1/2]. 

As in section |5j we fix 6, and fc* and express if? in terms of its near- and far-frequency compo- 
nents, around the band edge Eh, (k*): 

1p X = l/Wr + V>far = T _1 |V»near(fc)p6 t (x; fe) } + T _1 V'far.b (fe)p6 {x\ k)^ , (6.3) 
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where we define, for b = 0, 1, . 



Vw(fc) = X (I* - M < A r ) T fc . «> A }(fc) = x (I* - **l < A r ) (j*. (•, k),^(-,k) 



L 2 ([0,1]) 



L 2 ([0,1]) 

#ter,b(fc) = x(|fc-*.| > A r «5 6 ,, 6 )r 6t {V A }(fe) = x(l*-**l > \ r K,b)(p b {-,k)^{-,k) 
where 5i t j denotes Kronecker's delta function. Equivalently, one has 

r 1 / 2 /_ °° _ \ 

lp X (x) = / [ 1puea.T(k)Ub t {x; fc) + ^ 1ptar,b(k)Ub[x; fc) ] dfc. 

Recall that {pb(x; fc)}b>o form a complete orthonormal set in Lp Cr ([0, 1]), and satisfy 

(-(d x + 2mkf + Q(x))p b (x;k) = E b (k)p b (x;k), x g [0, l],p b (x + 1; k) = p b (x;k) . (6.4) 



Therefore, taking the inner product of (6.2| with p b (x;k), and using self-adjointness of the 
operator — (d x + 2-Kik) 2 + Q as well as the identity (6.4 1, yields 



(E b (k)-E x ) (p b (-,k),^(-,k)\ +x(p b (;k),(v^Y (•,*)) =0. (6.5) 

V / \ / L 2 ([0,1]) \ \ J I I,2([0,l]) 

or equivalently, using notation \2A\ , 

(E b {k) - £ A ) % |i/» A | (fc) + A7t, jV^} (fc) = 0. (6.6) 



We can now decompose equation ( 6.5 1 into near- and far-frequency equations, around E bt (fc*), 
the edge of the 6*-th band of the continuous spectrum. The coupled equations for T/> ne ar and i/>f ar 
read: 



[E b , (fc) - E x ) X {\k\< \ r ) (p b , (•,*), V A G, k)) 



1-2 ([0,1]) 



(6.7) 



+ A X (|fc| < A r ) (pi. (., fc), [V + vWP (•,*)) 



L2([0,l]) 



0. 



and for ieN: 



(£ 6 (fc)- J B A ) X (|fc| > A%„, 6 )^ 6 (-,fc), </>*(■, fc) 



£ 2 ([o,i]) 



(6.8) 



+ A X (|fc| > A r 5 6 ,, 6 )(p t ,(-,fc),[l/(^ nea r+^a 1 -F(-,fc)) i2([oa]) = 0. 

Equivalently, we write the near and far frequency equations in the form 

(E b ,{k) - £ A ) ^„ear(ft) + A X (|*| < A r ) (Tfe, {linear} (fc) + Tb, {V^} (fc)) = 0, 



(6.9) 



(E b {k) - £ A ) ^far, b (fc) + A X (|*| > \ r K,b) {Tb {^Vnear} (fc) + Tb {V^far} (fc)) = 0. (6.10) 

Equations (6.9 1 and ( 6.10 1 are, for the case of non-trivial periodic potentials, Q(x), the analogue 
of OH5.5L 
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6.2 Analysis of the far frequency Floquet-Bloch components 

In this section we study the far frequency equation ( |6.f0[ ). We will show that we can write it in 
terms of the near frequency solution and will determine a bound on the far solution in terms of 
the near solution. The next result is therefore the analogue of Proposition |5.1| and facilitates the 
reduction of the eigenvalue problem to a closed equation for the near-frequency components of the 
eigenstate. 

For clarity of presentation and without any loss of generality, we assume henceforth that we are 
localizing near the lowermost end point of the b*-th band and that fe» = 0. Thus, by Lemma \2.S\ 

&* is even, fc, = 0, with Eb„(0) = E*. 

N.B. For fc, — 0, note that Pb(x; fc,) = Ub(x; fc») and we use these expressions interchangeably. For 
fc* — 1/2 one has to distinguish between Pb(x; fc*) and Ub(x;k,). 

Proposition 6.1. Assume b t is even and consider E+ = E b , (0) the lowermost edge of the b t -th 
band. There exists Ao > 0, such that for < A < Ao, the following holds. Set 

>< A r - 1 i|^B i ,,(0)| 1/2 , 0<r<i. (6.11) 



Then, there is a unique solution ipf ar [linear, A] of the far-frequency system (6.101 in a small 
L 2 — neighborhood of 0. The mapping (^„„ r ; A) h-> Tpf ar = ipf ar [ipnear, A] maps L 2 (R) x (0,Ao) 
to H 2 and ipfar = 1 1 4>far satisfies the bound 

||y>/ar[l/Wr;A]|| H2(B) < C ( 1 1 V \ | i»o ) A 1 ~ 2r \ \ l/Wr || L 2 (R) . (6.12) 

Remark 6.2. Recall that we have assumed (1 + |a;|)l / (x) £ L X (R) and V £ L°° . It is in the 



proof of the bound (6.12 I that we have used V £ L°° . We believe it possible to work under milder 
assumption (1 + |a;|)l / (a;) £ L 1 (R). In this case, we would bound tpf ar in H 1 ^) and the analysis 
that would follow would be a bit more technical. We leave this an exercise. 

Proof. We begin by showing that there exists Ao > such that for all < A < Ao, there is a 
constant Ci > such that 

{Eb^-E^ > CiA 2r , A r < |fc| < 1/2 , (6.13) 
\E b {k) - £*| > Ci, b^b,, \k\ < 1/2 . (6.14) 



Note first that (6.141 is an immediate consequence of E* being the endpoint of the (6*) spectral 



gap. To prove ( 6.13 1 recall, by Lemma 2.2 that i5* = Eb, (0), an eigenvalue at the edge of a spectral 



gap, is simple, and k h-> Eb,_(k) — E+ is continuous. Therefore, for any Ao, such that < Ao < 1/2 

min \E b , (k) - EA> C*(A ) > 0. (6.15) 

A„<|fc|<l/2 ~ 

For |fc| < Ao, we approximate Eb, (k) by a Taylor expansion. In particular, since Eb, (fc) is smooth 
for fc near fc„ = 0, d k E b ,(0) = and dlE b ,{Q) / 0, we have E b ,{k) - E b ,(0) - \dlE b ,{Q) k 2 = 
C(|fc| 3 ). Therefore, we can choose Ao > sufficiently small so that for all A < Ao we have 

\E b ,(k) -E b , (0)| > ^ \dlE b ,(0)\ \ 2r , for all A < |fc| < Ao . (6.16) 



It follows from (6.151 and (6.161 that for sufficiently small Ao > 0, 



^>|fc|>A>0 ^ \E bt (k)-E*\ > mini ^\d'lE b , (0)| A 2r ,C(A ) 
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Thus if E x = E,- X 2 6 2 , 9 < A r_1 \ \8%E b , (0)| 1/2 , then for < A < A sufficiently small, there 
is a positive constant C\ , such that 



\E b ,(k)-E x \ > CiA 2r . 



By (6.131 and (6.14|, the far-frequency system, (6.101, may be re- written as 
x(l*l>A r * 6 . l6 ) 



tpiax,b(k) + A- 



G w [<t>; A] (*) 



E b (k) - E x 



T b {v^w} (*) + r b {vvfaj (*) ) = o, b > o. 

>L 2 by 



(6.17) 



(6.18) 



X (1*1 > A r (5b»,b) T b M (fc) + A x( ^y T t {V (0 + ^)} (fc) 



so that the triple (V'far, V'near; A) satisfies ( 6.18[ ) if and only if for all b G N, G(,^ nclI [V'fari A] = 0. 
Note that Gi,^ near [0; 0] = 0. Moreover, Gb,^ nesj . is smooth in a L 2 x R- neighborhood of (0;0). 
We will apply the Implicit Function Theorem to show that for all A in a small neighborhood of 0, 
there exists a unique solution = 0(A) such that Gf,,,/, noar [0; A] = 0. For this, it suffices to check 
that the differential, D^GfOjO], is an invertible map from L 2 to L 2 . We have 



X (1*1 >\ r 5 K , b )T b {/}(*) + A x g^j) 7-, {V (/ + ^ ncar)} (fc) 



[(^G 6 ,^ nrar [0;A])/] (*) = 
Thus, 

[(^G^ noar [0; 0]) /] (k) = % {/} (fc), 6 = 0, 1, 2, ... 
Now, {7& : L 2 ->• L°°([-l/2, l/2])} 6 N is invertible by the Floquet-Bloch decomposition: 

,1/2 



/(*) = £/ wx*)e 2 ™* 

i,_n -/-1/2 



Pb(x; fc)dfc. 



Therefore, there exists Ao > such that for all A < Ao, there exist a 0(A) G L 2 which solves 
Gb,^ noar [0(A); A] =0, < A < A . Denoting Vfar linear ; A] = 0(A), we obtain V6 G N, 
Gb,^ noar [V'farjA] = 0; thus (ipu T , line ar, A) satisfies ( |6.18 l. 

We now prove the bound (6.121 on ||V>far linear; A] ||h 2 (r)- Since ||fl||ji2 < C \\g\\ X 2 (Proposi- 
tion 



2.4|), it suffices to bound ||i/>f ar || A -2. First, we rewrite (6.181 as 

-A 



Tb {V'far [?/Wr; A]} = (1*1 > A^^b, ,b) Tb {V (tpneax + V^ar [V>near 5 A])} . 

Eb(K) — -£/ 



One has 



War [V'ncar; A] 



1/2 



^(l+6 2 ) 2 |T 6 {^ f ar}(fc)| 2 dfc 
2,2, 2X (1*1 > A r &.,») 



/1/2 
£(1 + & 2 ) 2 A 
-!/2 b= o 



near + Vfar [V'near! A])}| dfc. 
|£b(fc) - i? 
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Now, by ( pTft , for |fc| > A r one has \E b ,(k) - E x \ 
asymptotics to write | (E b (k) - E^ 1 1 ~ \{b 2 - E,)' 1 

.12 . r 1 / 2 

[V'far linear; A] 



< CiA~ 2r . For b b„, we use Weyl 
(b 2 + We therefore have 



/1/2 ^ 
V|Tfe{y(VVear+ V'far [Vn 
-V2 6=0 



;A])}| 2 dfc 



+ A 2 



1/2 



1/2 



(1 + 6, 2 ) 2 X (|*| > A'') \Tb, {V (Vncar + Vfar [Vnear ! A])}| 2 dfc 



< A 2 4r (||(V'V , near)~|| A rO + || (V^far [Vncar! A])~|| A . ) , 

Using Proposition |2.4| to return to H" norms, we have 

H^far [V'ncar; A] || H 2 < A 1 "^ ( 1 1 VVncar 1 1 L 2 + ||VV&r [V'ncar J A] || £2 ) 

< C(||V||loo) A 1 " 2 " (||Vncar|| L 2 + ||Vfar ^near! A]|| ff2 ) . 

Recall that < r < 1. If we further constrain r so that < r < 1/2, by choosing A small enough 
so that C(||V||z,«.) A 1 - 2 '' < 1/2, we obtain ||Vfar [Vw; A] || H2 < C (|| V||i,°c ) A 1 ~ 2r ||Vw|| L 2, for 
< r < 1/2. The proof of Proposition 6.1 is complete. □ 

6.3 Analysis of the near frequency Floquet-Bloch component 



With the properties of the map tp n 



Vfar [Vnear, A] now understood via Proposition RTT] we now 



view and study (6.9 1 as a closed eigenvalue problem for (E , ip n 



E b . (fc) - E*j Vncar(fc) + \ Xx r (fc)7S, { V^ncar} (*) + A X X r {k)Tb, {^Vfar [V'ncar; A] } (fc) = 0. 

(6.19) 

Equation (6.191 is localized in the region \k\ < A', < r < 1/2. By careful expansion and 
rescaling of (6.191 we shall obtain an equation, which at leading order in A, is a perturbation 



of the general class of equations to which Lemma |4.1| applies. The size of the perturbation is 
estimated in Proposition |6.6| and the perturbed equation is then solved by applying Lemma |4.3| 



Theorem, making explicit the leading and higher order contributions. 



In Lemmata 6.3 6.4 and |6.5| we expand the three terms in (6.191 about fc* = using Taylor's 

There exists k' 



Lemma 6.3. Denote E x 
such that \k'\ < A r , and 



E, - A 



Eb„ (0) — X 6 , as in Proposition 



6.1 



E b ,{k)-E 



v )Vw(£;) = Qdfc-Eb.tO) k 2 + X 2 6 2 ^j Vw(fc) + A-Ro [Vw;A] (*,*'), 



where 



■0 [^near! A] (k,k') = j ^ fc 4 8^E K {k') Vncar(fc). (6.21)) 



Proof. Taylor expansion of E bt (fc) about fc* = to fourth order and use of E x = E b , (0) — A 2 # 2 
and d J k E b ,(0) = for j = 1, 3 we obtain E b ,{k) - E x = \d 2 k E b ,{Q)k 2 + X 2 6 2 + ^diE b , (fc')fc 4 , 
which is equivalent to (6.201. □ 



Lemma 6.4. There exists k" such that \k"\ < X r , and 

%, {VlP n ear}(k) = (pb, (■ | 0) , p b , (■ | 0)T { V T~ 1 {Vncar}} (-, *) 



i 2 ([0,l]) 



+ i?l 



V'ncar; A J (fc, fc"), 



2G 



with 

Ri [</w;A] (M") = <Pfc,0;O) I T{yfi}0 J ft)> i2([o , 1]) 



+ fc (d ft p 6 , (•; k"), Pb , (•; 0)T {V-F^^ear}} (•. *)) ta + fc <^Pi>, hO-T^} (•, fc )) i 2 {[0 , 1] 

(6.21) 

iti/iere £1 = T -1 jv>near(fc) (pb„ (a;; fc) - Pi,, (a;; 0)) } . 



Proof. Let us recall that by definition (6.3 1, one has 



1pnezr(x) = j ^near (• )p&, (»; ■)} • (6.22) 

Since V'near(fc) = x(\M < A r ) V'near(fc), we decompose: 

^near(a:) = {V'nearOPb. 0^ ')} 0*0 = Pf>, °) F~ 1 { V'near } + £l {x) (6.23) 

where 

= T" 1 {^ncar(-)(P6,(a;;-) ~ Pi>, 0)) } ■ (6.24) 

Above, we used that T -1 commutes with multiplication by a 1- periodic function of x, and that 
when acting on a function which is localized near k = 0, and which does not depend on x, T _1 is 
equivalent to the standard inverse Fourier transform; see secti on [2] 

The decomposition of Tb, {VV>near} (k) follows now from (6.231 and the Taylor expansion of 
Pb, (x, k) around fc, = 0: 

p b , (x; k) = p b , (x; 0) + kd k p b , (x, k"), (6.25) 
where \k"\ < X r , depends on k. 

%, {Vlpnear} (fc) = (Pb, (• ! k) , T { Vl/Wr } (• , fc)) L 2 ([0jl]) 

near} fc)) £, 2 ([0,1]) + (kdkPb, (s fc ) , T { Vl/'near } (■ , fc)) Z ,2Q 01 ]j 
= (p bt (-;0),p i) ,(-;0)r{VJ- 1 {^nea r }}(-;fc)) L2([oi]) + «1 [^near! A] (fc, fc")- 

Here, R\ ^ nC ar; A J is defined by: 

Rl [V'ncar; Aj (fc, fc") = (p&, (•; 0) , T {V£ 1 } (■,fc)) i 2 ([0 , 1]) +fc(a fc p i) ,(-;fc"),T{Vl/' ne a r } (') *)) i 3([ 0|1] ) , 

and one obtains ( |6.21[ ), using once again the decomposition ( |6.23[ ) on the second term. The proof 
of Lemma |6.4| is now complete. □ 

We next give a precise expansion of the leading order term in Lemma |6.4| 
Lemma 6.5. For any \k\ < X r , there exists fc'" with |fc"'| < 2A r , such that 



{p b , (•; o), j*. (•; o)T {kf-^Vw}} (•; *)) (6.26) 

|p b , (a;;0)| a V(a;)da:J / rp near (l)dl + R 2 ^neor] {k,k"'), 



with 
Rz 



linear] (k,k"') = -2m f \p b , (x; 0)\ 2 xV (x)e- 27rik "' x dx [ (k - l)i> near (l) dl. (6.27) 

J — oo -/ — oo 
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Proof. By the definition of T, one has 

T { KF" 1 Ww}} (x; k) = Y l e 2nlnx T { V-F-^Vw}} (* + n) = ^ e 2 ™* V(fe + n - 1)&«,(Q dl. 



n6Z nGZ 

Therefore, one has the identity 

(pb,{-,0), Pb ,(-,0)T{vF- 1 {^ ncal }\(-,k)) = Vp(n) / ^(fc + n-J)Vw(0^, 

\ L J /^Oo.H) t& J -°° 

where we defined T{n) = (p b , (■, 0) jP{ ,, (-, 0)e 2 ™-) i2([01]) . 

Since |fc| < A r and |/| < A r , one has \k - l\ < \k\ + \l\ < 2\ r . We Taylor expand V(k + n - I) 
as follows: 

V(k + n - I) = V(n) + (k- 1)9' (k'" + n), 
for some k'" such that \k"'\ < 2A ? . The first term of (6.261 now follows from the identity: 

V(n)V(n) = (po* (S 0),p b * (■; °) e2 ™') L2 . I0 m V(n) = f W {x; 0)| 2 £ e™*V(n)dx 

roc 

\p b *(x;0)\ 2 V(x + n)dx = / \p b * (x; 0)\ 2 V(x)dx. 



Here, we used the Poisson summation formula and that x h-> p;,* (a:; 0) is 1— periodic. 
Similarly, once we note that 

V'(k"' + n) = F[ - 2-Ki{-)V{-)e- 2lTlk "'-} (n), 

we can show as above that 

-2-Kik'" x 



/oc 
\ Pb *(x;0)\ 2 xV(x)e-^' K x dr. 
-oo 



This completes the proof of Lemma |6.5| □ 

The rescaled closed equation. Using Lemmata |6.3| |6.4| and |6.5| one can express the near 
frequency equation (6.191 as follows: 

(^£ 6 ,(0)fc 2 + AV)Vw r (fc) + A X (|fc|<A r )(y \ Pb ,{x;0)\ 2 V(x)dx^ J ^(l)dl 

= -A X (|*| < A r )7l[Vw;A] (k), (6.28) 

where TZ [i/w; A] (k) = T b , {Vipt^} + R + Ri + R 2 - 

Seeking to extract the dominant and higher order terms in A, we introduce the scaled near- 
frequency components: 

&•*•(*) = i$A f |j = («) , where k = A k . (6.29) 



Expressing (6.281 in terms of <E>a and k we obtain, after dividing out by A, 

= -x (|/s| < A''" 1 ) n [Vw; A] (A/s) = i?($ A ). (6.30) 
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Equation ( 6.30 1 is of the form £.[8]&\(k) = _R($a), where C[6] is given by ( 4.1 1 with parameters 



A= ^dtE b ,(Q) , B = - I \p b ,(x;0)\ 2 V(x)dx , and = 1 - r 



In order to solve (6.301 via Lemma 4.3 we need a bound on _R(«I>a) of the form (4.191. 

Proposition 6.6. Assume that V is such that (1 + | ■ QV(-) G L 1 and V € L°° . Then R(3>\), 
defined in (6.301, satisfies the bound 



R($x)\\ l2 i = \\x (I • I < A*- 1 ) Tl [i> near ; A] (AOlLa.-x < C\ 



a(r) 



1 1,2,1 



(6.31) 



where a(r) — max{| — 2r,2r, }• The cons t an t, C depends on: 



Mlx,i'IML°°' sup ^|(j3 6 »(-;fe), 



'£ 2 ([0,1]) I 



sup I d%Eb,(k)\, sup ||9fePk,,(-;fc)|| J 



|fc|<A>- 



|fe|<A r 



\pb, (x; Q)\\V(x)\ dx and i \pt,*(x\0)\ |a;V(a;)| dx, and is finite by Lemmata 



2.3 and 



2.5 



Proof of ProposiUon \6.6\ Recall that TZ(Xk), the right hand side of (6.301 has the form 

72. [iptu linear; A] , ^> near ; A] (A«) = x < A r_1 ) r7fc, { V^f ar } (A/t) + i? [fcar; aJ (Ak, fc') 

+ i?l ^near; A| (A«, fc") + i?2 [^nearl (Aft, fc'")) 

= (/) + (II) + (III) + (IV) . (6.32) 

We proceed by estimating each of the terms: (J), (II), (III) and (IV). 
(I) Estimation of X (M < A r_1 ) T b , {Vt/w} (Ak): We have 

x(\ K \ < A r_1 ) 



|x(l-l<A I - 1 )r 6t {^ far }(A.)|r L2 ,„ 1 = 



1 + K 2 



< \\Tb, {V^r}\[ 



\Tb, {FVfar}(AK)| 2 dK 



We now consider 71, {Vipfa,r} (■) in detail. By definition, one has 
%, {V^}(k) = { Pb ,(-,k),T{V^}(;k)) LH[0A]) 

/oo 
V(k + n-l)^(l)dl 
- oo 



= { Pb 



^2\Pb,(-;k),e 2 



V(k + n-l) i 
LHlo,i])J- x (l + WY' 2 ' 



L 2 ([0,1]) 



(l + \l\ 2 ) 1/2 ^ r (l)dl. 



Moreover, 



f°° V(k- 



+ n-l) 



|/| 2 )!/2 



(l + \l\ 2 ) 1/2 4> iai (l)dl\ < \\v\\ TOO HVfatll^ < A^IVwII 



< A 



I £2 



A 1_2r A~5 S A , 



II 2 



where we used Proposition 6.1 definition (6.291 and, by Proposition 

|2 



2.4 



^near 1 1 £2 



T 1 {V'ncar(fc)p6»(x; fc)}|^ 2 < j | l/> nca r (fc)p 6 , (iE; fc) | 



A"' 



/ j^noar(A:)| dfc 



-1/2 



(6.33) 
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Finally, it follows 



\\T h . {V^ r }\\ LX , < A5- 2r C 11$ 



(6.34) 



with C = C ^\V\\ LOO , \\V\\ LOO , sup |fc|<A . £ nez I (Pt, (■; k), e 2 ™'> i2([0il]) |) ■ 

(77) Estimation of \ < A'" -1 ) i?o jV'nearjAj (Xk, k'), give in (6.201: We have (constants 
plicit) 



Z I — K I ^ I a 

= A 4 / X (|«|<A r - 1 ) * A ( K ) dK 



Pi Z aW X(l«| < A r - 1 )(l + K 2 ) S A («) d«<A 4 sup 

-00 I 1 + K J 1 1 |«|<A'— 1 



(1 + K 2 



1 1,2,1 



< A 4 



Therefore, 



|x(N < A r 4 ) 7?Q [V'ncarjAj (A/C, fc') j 



£2.-1 



x(l«l < A r_1 ) ia£Efc.(*')*V*A(«) 



< A 2r sup |to.(fc')| ||Sa| 



\k'\<\ r 



£2.1 



£2.-1 

(6.35) 



(777) Estimation of X (|«| < A r_1 ) 7?i [v 
Hi 



A (Ak), given in (6.211: Recall 



/£ 2 ([0,1]J 



Vnear; a] (fc, *") = ( Pb , {-, 0),T{V£ 1 } (■, fc)>, 

+ fc (■; k"), Pb . (•; 0)T {V^'M^ear}) (•, k)\ 

\ 1)1 £ 2 ([0,1]) 



+ k(d k p b ,(-;k"),T{V£ 1 }(;k)) 1 



^£2([0,l]) ' 

where B\ = T" 1 |V" nC a r (fc) (x; ft) - p 6 , (as; 0))|. 

Let us first obtain an estimate on £\. Using Taylor expansion (6.251, one has 

fl/2 



(6.36) 



\£,{x)\ = 



r 1 / 2 _ 

/ e 2mkx tp neiiI (k) (p b „ (x; k) - p b , (x; 0)) dk 

J -1/2 



< sup \dkPb,(x;k' : 

x£M,\k"\<\ T 



\kx(\k\ < X r ) ?/>near(fe)| dk 



< X sup \\dkPb,(-;k")\\Loc / |kx(N<A'' 4 ) $a(^)| 

fc"|<A'" 7-oo 



< A sup \\d k pb,(-;k")\\ L ^ ( / — 

\k"\<X T \J\k.\<\t-~l 1 + 



dK 



1/2 



I £3.1 



< 2A^ sup ||«kP6.(-;*")lk- 

|fc"|<A»- 11 



£2.1 



so that we deduce 



\\£x\\ LOO < 2X^ sup \\d kPb ,(x;k")\\ L ~ 11$, 



*"|<A r 



I £2,1 



(6.37) 
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Estimation of the first term of (6.361 is as follows. One has 



\x{\k\ < A r - 1 )(p fc .(.;0) i r{V£i}(. ) A«)> ia([0il]) |^_ 1 = 

X(I ^ A 2 - \{ Pb A-,0),T{V£ 1 }(;\K)) LH[0A]) \ 2 dK. 

Turning to the integrand of the above expression, we rewrite the inner product 



<Pb.(sO),T{V£i}(-,A/s)> ia([0 , 1]) = / T{p b ,(-,0)£i(-M-)}(x;\ K ) 

Jo 

= / X) e2nlnXjr {P^ (■> 0)Si(-)V(-)} (A* + ") 

= J 7 {p».(-;O)£i(-)V(0}(AK), 
where we used that pt, (x; 0) is 1-periodic, so that it commutes with T, and the Poisson sum- 
mation formula. It follows that \{pb t (•; 0), T{V£i} (•, Ak))^^ ^ J < ||pi„ (■; 0)£i(-)V r (-)|| Ll < 
|£i ioo / \Pb* { x \ 0)||V(x)| Using (6.371, one deduces 



|x(l«l<A r - 1 )<p*.(-;0),r{V r e 1 }(.,A«)> ia([Oil] )|| ri?i _ 1 < CA^|$ a ( k )| l2i , (6.38) 



with C = C(su P|fe „ |<A . \\d kPb , (■; fc")|U- , / K (at; 0)||^(x)| dx). 
We can bound the inner product of the second term of (6.361 by 



d kPb , (■ ; k") , Pi. (■ ; 0)T ( U^" 1 { Vncar U (•, A«)\ 

I- I J J / £2([0,1]) 

9fcP6» (x; ( X ; 0) ^ g 2 ™^ J r {VJ'" 1 {^near}}(AK + Tl)dx 

/ d k pb,{x\k'')p bt (x;0)V(x + n)J r ~ 1 {^ nca , r }(x + n)e~ 2nl<XK+n)x 
nez Jo 

/OO 
9fcPt, (x; k")p b , (x; fyVix)?^ 1 {ii ne ^}(x)e~~ 27TlX, * x dx 
- OO 

/oo 
(a^;fc )pt* (a;; 0)V (x) \dx 
-oo 



ci.i.- 



< Vnear Sup \\d k p b ,(-\k")\\ L o 

II II i,i |V|<A r 



|pb,(a;;0)| |V(x)|da;, 



where we used the Poisson summation formula and the periodicity properties of p br and d k p bt 
Now, note that 



so that 



/oo poo P OO -| 1 1 1 1 

\^r(l)\dl = J \$x(v)\ dr, = y - (1 + r? 2 ) 1/2 j$ A (77)1 rfq < C 

X (|/S| < A r_1 ) Aft (&Pi. (•; fe"),Pi» (•; 0)T { V^H^near}} (•, A«)) 

< ca||$a(k 



IL2.1 



L2([0,l]) 

« 2 X (M < X^ 1 ) \ 1+r ||^ , 

U ' 1] < AT- * A (/s) 



1 +K 2 



(6.39) 



31 



with C = C (sup |fc „ |<Ar \\d k p b „(-,k")\\ L °o, J\ Pbt (x;0)\\V(x)\dxy 
The last term in (6.361 is estimated similarly. Note that 

|(<9 fc p { ,, > (-;fc"),T{V£i} (-, Ak)) l2{[01]) | 

= / dkPb, (x; k") e 27V ' nx J r {V£i}(\K + n)da 

/OO 
dkPb, (x; k")V(x)£i(x)e~ 2lTt KX dx 
- OO 

< ||£i||z,«> ||V|| £ i sup ||9fcP6»(s*")|U°° 

|k"|<V 

\ 2 

ll#A 



< A— ||V|| £1 ( sup \\d kPb ,(-,k")\\ L o 

|k"|<A*" 



where we used the Poisson summation formula one more time along with the periodicity of dtpt, 
and the bound ( |6.37| |. Therefore, one obtains 



|x (N < a*- 1 ) ak (a fcPb , (■; fc"),r{vfi} (., A«)) i2([0] 

where C = C (\\V\\ L i, sup !fc „| <Ar \\d k pb, (•; &")IU°°) • 
Finally, we have 



<CA i+r 



(6.40) 



X(|«| < A r_1 )iii U near ;A (Ak) < CAT $J 



(6.41) 



withC = C ||V|| £ i, sup ||flkP5.(-;*")L-. / \Pb*&0)\\V(x)\ dx 

\ |fc"|<A r J 

(IV) Estimation of X (M < A'' -1 ) i? 2 [^near](AK, fe'"), given in ( |6~27| ): Recall 



fl2[^neat](AK,fc'") = T(k) x 2th / |p 6 . (z; 0)\ 2 xV{x)e 



dx 



where 



X( K ) = -x(M < A 1 - 1 ) / (A(«-^))x(M<A p - 1 )* A (T 7 )di/ 

•/ — OO 

The integral, is bounded in L 2,_1 (1R) as follows: 



(6.42) 



imi 



£2,-1 



< A 



1 + K 2 

x(H <A r - x ) 

1 + K 2 



(k - 77)x (M < A r x ) $Ad?7 



A 2 $A 



M<A r 

(^-r,) 2 



1 + 77 S 



-drjdn KE>a 



£2,1 



^^J^Jr, (1+K 2 )(1 + ?? 2 ) 



X (|«|<A r - 1 )x(W<A'- 1 )d«d 7? . 



One easily checks that 



Z / (i + « : 



-?7) 2 

)(l + r, 2 ) 



x(|«|<A r - 1 )x(W<A r - 1 )d«^ < X r 
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so that 



x(|«| < A r - 1 ) J R 2 [Vnca r ](A«,fc'") a _ 1 <C / \ Pb *(x;0)\ 2 \xV(x)\dx) U x \\ 2i A~ 



Altogether, (6.341, (6.351, (6.41 1, and (6.431 yield the estimate ol Proposition 6.6 



(6.43) 

□ 



6.4 Completion of the proof of Theorem 3.4 



We now prove Theorem [3^4] by an application of Lemma |4.3| to equation ( |6.30[ ), where the remainder 
is estimated in Proposition |6.6| 

Proof of Theorem\3^ We seek E x = E bt (0) - X 2 9 2 and ip x of the form 

( OC 

1p X = ^noar + Ipiat = T~ 1 j V'neaf (fc)Pf>, (x;k) j + T~ 1 < ^ ^far ,b (fc)?>ft (zj fc) 
1/2 
-1/2 



where i/'near, V'far satisfy equations ( 6.9 I— ( 6.10 1 ; see section 6.1 By application of Proposition 6.1 
one has that ^f ar is uniquely defined as a function of y> nC ar and A, and that H^tar linear; A]|| h2 < 
A 1 "" 



-'near \ \ 



Then, defining $ A as in (6.291: 



Vw(fc) = j®\ ( j ) = t $a(k) 



1 



A ' 



(6.44) 



By Proposition 6.6 the rescaled (from (6.9l) near-frequency equation (6.301 can be written as 



8%E b .(0)K 2 +6 2 X xr -i(«)*A(«) + X xr -i(«) / K(-;0)|V / JC^fa)^)*? 



"X (N < A r_1 ) TZ ($a) (^), (6.45) 



i-+i ■ 



with ||"ft($ A )|| i2 ,-i < CX a(r) ||$ A || J . JI>I) anda(r)=max(i-2r,2r, 

From now on, we set r = 1/8, a = 1/4, which yield optimal estimates. Applying Lemma 4.3 
with /3 = 1 - r = 7/8, 



.4 



8tt 2 



dlE b ,(0) and B 



(x; 0)\ 2 V(x)dx ^assumed to be positive^ , (6.46) 



we deduce that there exists a solution ^ 2 , <1>a ) of the rescaled near-frequency equation (6.451, 
satisfying 

ll*A-/o|L-[ < C* AJ and \8 2 - 8 2 \ < C A' . (6.47) 
Here ^#o(A),/oj is a solution of the homogeneous equation 

£o,a (0o,fo) =(4n 2 Ae+e 2 )f ~B X (\!i\<X-i) f°° % (\v\ < To(v)dr] = 0, 
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as described in Lemma 4.1 



Thus vWr(£) = (|) and E x = E b , (O)-A 2 2 (A) are well-defined 

(and satisfy the Ansatz of Lemma 6.31, and i/>f ar is uniquely determined as the solution of ( 6.10 1; 
see Lemma [6JJ It follows that 



lp X (x) = ^far + tp n 



V'far + 



1/2 



VWr(fc) M S<» (x; k) dk 



(6.48) 



1/2 



is well-defined. 

There remains to prove estimates ( [3~9| and ( |3.10[ ). Recalling that E x = E b , (0) - A 2 6> 2 , | |6.47 l 
implies £ A - (E b ,(0) - A 2 6»g) | < CA 2+1/4 . By Lemma [Zlj one has j(9 (A) - < 



C(A, _B)As , so that one can set 



-B 



-\fZ \^(pr,k.)\ 2 V(x)d a 



4/1 



£*dlE K (k* 



and estimate (3.9 1 follows. 



We now turn to a proof of the eigenfunction approximation (3.101. Recall 

"1/2 _ 



ipne^(k)u br (x;k) = 



-1/2 
1/2A 



1/2 



1/2 



A ' 



(x; fc) dk 



1/2A 



x(lCI < A-5)$ A (|)e 2 ™*%,,(z;A6d| 
j X (ill < A"') $a (|) e^-^. (as; 0) ^ 

+ / X (l€l < A"*) $ A (Oe^^AO^^fc'X 



«6»(x;0) / x III < A" 5 /o(|)e : 



d| 



+ «b.(s;0) [ x(\Z\<^ i ) (*a-/o) (Ik 



+ / x(|ll<A-i)$A(|)e 27riA£a: (Ae)5 fc p 6 ,( a; ;fc')^ 
= /i(*) + I 2 (x) + /3(a;) 1 
with |fc'| = |fe'(A|)| < As. Now, since X (ill < A~s) /o(£) = Ml), one has 

h(x) = u^ix-O)^ 1 {£>} (Xx). 



By (4.81 in Lemma 4.1 one has 
2 



sup 



h{x) - — u b , (x; 0) exp 



AS 



i l6 ,(x;0)|^ 1 {/ o }(A a; ) 
< C|K(-;0)]|^ A 7/8 . 



sup 



exp 



as, : 

2A N 
(6.49) 



and ||pb» (•; 0) is bounded; see Lemma 2.5 
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Let us now estimate Iz(x) and Is(x). One has 



\h(x)\ 



u b , (x; 0) J X < A"*) (* A - /o) (Oe^' * 
r xfl£l < A - ^ , 

< lPb.(»;0)l / , , m , / / (i + |{| 3 ) 1/ 1^)-/o(fl 



< C\\p k (-,0)\\l°° ||$a-/o|| < C(yl,B)||p i) ,(-;0)||i,-A 1/4 , 



(6.50) 



where the last inequality comes from (6.471. Similarly, 

|/3(2)l = ^x(lCI< A-s)$ A ( ? ) e 2 ^(AOa fe p 6t (x;fc') ^ 

<A sup ||ftkPfc.(-;*')lli» /"x(|CI<A~J)|€||$a(0|* 



< C sup \\d k Pb,{-\k')\\ L o 
|k'|<AVa 



Hi? 



(6.51) 



By ( |6.49[ ), (|6_50j and ([6_51], one has 

2 



/i(ar)+/a(i)+/ 3 (a!) 



B Ui,„(a;;0)exp ^ 



AS, 



x\ ) + V<rem(a;), with HV'remlU^ < A 1/4 . 



Finally, let us note that by Sobolev embeddings, one has 

||iMz,« < IIVtar||fl» < ca 1 - 174 !!^!!^ =ca 1/2 - 1/4 ||*a|U2 <ca 1/4 , 



where we use Proposition 6.1 with r = 1/8, and (6.331 



It follows that ip = ipncar + ^far satisfies 
2 



sup 

16M 



ip x (x) 



B 



Ub, {x; 0) exp(Aao|a;|) 



< CA 1//4 , with ao 



B_ 
2A' 



Since ip x is defined up to a multiplicative constant, (3.101 holds. This completes the proof of 
Theorem ISH □ 



A General properties of E^k) and derivatives c^-E^fe*), 
where Eb(k*) is the endpoint of a spectral band 



For the sake of completeness, we prove Lemma [2. 2| which concerns the spectrum of the eigenvalue 
problem, for E fixed, 



{-dl + Q{xj) iP(x; E) = Ej>(x; E), Q(x + 1) = Q(x), 
with solutions which satisfy 

i[>(x + 1; E) = pip(x; E) peC. 



(A.l) 



Let 4>i(x; E) and <j>2(x; E) be two linearly independent solutions of (A.l \ such that 



M0;E) = 1, 



<t>' 2 (0;E) = l. 
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The functions fa(x + 1; E) and 4>2(x + 1; E) are two other linearly independent solutions to ( A.l I, 
so that we can write 

fa (x + 1;E) = A n fa(x; E) + A 12 fa(x; E), (A.2) 
fa(x + 1;E) = A 21 fa(x; E) + A 22 fa(x; E). (A.3) 



Note that the matrix (Aij) is nonsingular. In general, every solution of (A.l I has the form 

i>(x; E) = cifa(x; E) + c 2 fa(x- E). (A.4) 

As we are specifically interested in solutions which satisfy if)(x + 1; E) — pip(x;E), one has the 
following identity 

ip(x + 1;E) = pfax- E) ^ a(fa(x + l;E)- pfa(x; E)) + c 2 (fa{x + l;E)- pfa(x; E)) = 

<^(d (An -p) + c 2 A 2 i) fa (x; E) + (ciA 12 + c 2 (A 22 - p)) 4> 2 (x; E) = 

' ci (An - p) + c 2 A 21 = 
ciAw + c 2 (A 22 - p) = 



The solvability condition ( |A.5[ ) is satisfied for nontrivial ci and c 2 if 

det(A - pJ) = 0, i.e. p 2 - (An + A 22 ) p + det(A) = 0. (A.6) 

Using that the Wronskian, W [<j>i , <f>2\ (x; E) = <f>i(x; E)(f>' 2 (x\ E) — 4>'i(x; E)(f>i(x; E), is constant 
with respect to x, one has 

det(A) = W l<t>i,H (1; E) = W [fa, cj> 2 ] (0; E) = 1. 
Therefore p must satisfy p 2 — D(E)p +1 = 0, where we define the discriminant 

D{E) = A 11 +A 22 = fa(l;E) + <j>' 2 (l;E) . (A.7) 

We note that the two solutions of the equation p 2 — D(E)p +1 = satisfy jpj < 1 if and only 
if the discriminant \D(E)\ < 2. In that case, one can write p = e ±2,rik , with k € (-1/2, 1/2], and 

D(E) = 2cos(27rfc). (A.8) 



As \p\ = 1, ip(x\ E) is a bounded solution to (A.l I, and E = Et(k) is in the continuous spectrum 
of Hq = —-^ri + Q- More precisely, for E — Eb(k), one has 

ijj(x;E h (k)) = u b (x;k) = e 2 " x p b (x;k), p b {x + 1; k) = p b (x; k). 

where {E b (k),p b (x; k)} b>0 , is the eigenpair solution to ( |2.2[ ), as defined in Section |5J For each 
6GN the eigenvalues {E b (k), k e (—1/2, 1/2]} sweep out the b th stability band B b = [G b> F b ]; see 
figure [2] 

Let us now rewrite Lemma [2,2| which states some of the properties associated with the stability 
bands. 



Lemma A.l (Lemma 2. 2 1. Assume E b (k*) is an endpoint of a spectral band of —d 2 + Q(x), which 
borders on a spectral gap. Then k* € {0, 1/2} and the following results hold: 

1. b even: E b (0) corresponds to the left (lowermost) endpoint of the band, 

E b (l/2) corresponds to the right (uppermost) end point, 
b odd: E b (Q) corresponds to the right (uppermost) end point of the band, 
E b (l/2) corresponds to the left (lowermost) end point. 



3G 




Figure 2: Discriminant, D(E). 



dkEbik*) = 0; 

b even: d 2 k E b (Q) > 0, 9^(1/2) < 0; 

b odd: d 2 k E b {Q) < 0, d 2 k E h (l/2) > 0; 
d s k E b {k,) = 0; 

E b (k t ) is a simple eigenvalue of the eigenvalue problem (2.1 1. 



The proof of Lemma A. 1 is a consequence of the following result, concerning the problem ( A.l I, 
and which is proved in the first two chapters of [101 and part I of |21j . 



Theorem A. 2. Consider the equation (A.l I, and define D(E) with (A.7l. Denote the edges of 
the stability bands as 

Go < F < F 1 < Gi < G 2 < F 2 < F 3 < G 3 . . . 
Then the following facts hold: 

I In the interval [G<irn, Fam], F>(E) decreases from 2 to —2. 
f In the interval {G 2 m,F 2m ), D'(E) < 0. 

II In the interval [i^m+i, Gim+i], D(E) increases from —2 to 2. 
II' In the interval (F 2m +i,G2 m +i), D'(E) > 0. 

III In (-oo,G ) and (G 2m+ i,G 2m+2 ), D(E) > 2. 

IV In (F 2m ,F 2m+1 ), D(F) < -2. 

V D(E) = ±2 and D'(E) = if and only if E is a double eigenvalue. Furthermore, D"(E) < 
ifD(E) = 2 and D"{E) > if D(E) = -2. 

Proof of Lemma A.l Let us recall that one has from (A.8l that the discriminant satisfies 



D(Eb(k)) — 2cos(27rfc). It follows that as k increases continuously from to 1/2, D{E) decreases 
continuously from 2 to -2. Therefore by I and II, G 2m (k) increases continuously from G 2m to F 2rn 
while G 2 m+i(k) decreases continuously from G 2m +i to F 2m +i. This proves claim 1. 

We skip to part 5 to show that if Eb(k,) corresponds to a band edge, that is if there exists a gap 
between the 6 th band and the closest consecutive one, then Eb(k„) is a simple eigenvalue. Without 
loss of generality, we assume Eb(k») to be the lowermost edge of an even band, for example G 2m 
in figure [2] Therefore, for any 5 > sufficiently small, 



D(Eb(k«) -S)>2 and D(E b (k*) + S) < 2. 



(A.9) 
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Assume for the sake of contradiction that Eb(k*) is a double eigenvalue, which means, by part V 
of Theorem A.2 that D'{E b {k t )) = and D" '(E b (k*)) < 0. Now, Taylor expand the discriminant 
about E b (k*), 



D{E) = D(E b (k*)) + D'(E b (k*))(E - E b (k*)) + -D" (E b (k t ))(E - E b (k,)) 2 +0((E- E b (k,)f) 
= -2 + ^D"(E b (k,))(E - E b (k,)) 2 +0((E- E b (K)) 3 ) . 

Since D"(E b (k*)) < 0, we have D(E b (k,) - 5) w 2 + (l/2)D"(E b (k,))8 2 < 2, which is a contra- 
diction of (A.9l. Therefore part 5 is proven and we have that at the band edges, E b (k*), the 
derivative of the discriminant is nonzero, 



dD 
dE 



(E b (k t )) ± 0. 



(A.10) 



To see part 2, note that — 47rsin(27rfc*) = dkD — 8eD dkE b . Using (A.lOl, we conclude that 
d k E b (k) = if and only if k = or k = 1/2. 

To prove part 3, we take the second derivative of D(E b ) with respect to k and evaluate at fc„, 



-87r cos(27rfc*) 



d 2 D 

~dk? 



(fc. 



d 2 D (dE b 
dEj \~dk 



(fc.) + 



dE b dk 2 



dE b dk 2 



•(*.) 



Therefore, by I', II', and ((ATO} we conclude 3. 

Similarly, to show 4, we take the third order derivative of D(E b ) with respect to k, 



d?D 



?6 d 2 DdE b d 2 E b 



dE b dk 3 



+ 3 



d 3 D /dfii 



dE 2 dk dk 2 



dE 3 



dk 



Evaluated at we have = ^S'i^*) ~ if§~ d dk^ w hich concludes the proof of Lemma 
once we again note (A. 10 1. 



2.2 

□ 



B Regularity of k i— > ^(fc) and /c i— >• /c) 

In this section we give a self-contained discussion the of regularity with respect to k of the Floquet- 
Bloch eigenvalues and eigenstates. 

Consider the k— pseudo-periodic eigenvalue problem for each k £ (—1/2, 1/2]: 

(— d 2 + Q(x)) u(x; k) = Eu(x;k), u{x + 1; k) = e 27rlk u(x; k) (B.l) 

Introducing the Floquet-Bloch phase explicitly via u(x; k) — e 2ntkx p(x] k), we obtain the equiv- 
alent formulation 

H Q {k)p(x;k) = (-(d x + 2Tvik) 2 + Q(x))p(x;k) = Ep{x;k), p{x + 1; k) = p(x; k) . (B.2) 



As noted, for each k £ (—1/2, 1/2], the eigenvalue problem ( B.2 1 (equivalently (B.l |) has a discrete 
sequence of eigenvalues Eo(k) < E\{k) < E2(k) < • • • < E n (k) < • • • . 

It can be proved, using the min-max characterization of eigenvalues of a self-adjoint operator 
that the maps k h-> E b (k), b — 0,1, ... , are locally Lipschitz continuous. A proof based on standard 
perturbation follows from results in |25j . An elementary proof is given in Appendix A of [llj . 

In the present paper, we require a Taylor expansion of the E b (k) near k = k», for which E b (k*) 
is the endpoint of a spectral band, which borders on a spectral gap. By part 5 of Lemma [A. 1[ the 
eigenvalue E b (k t ) is simple. We prove the following 
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Theorem B.l. Suppose is the endpoint of a spectral band of —d 2 + Q(x), which borders on a 
gap. Thus, = E b (k t ) for fc* € {0,1/2} and the corresponding eigenspace of solutions to (B.2I 
has dimension equal to 1. We denote the normalized eigenfunction by p(x;k„); 

\p(y;k,)\ 2 dy = 1. 

Then, there exists p > such that for all complex k in a complex disc centered at fc*, B p (k t ) = 
{k G C : |fc — As«| < p}, the following holds: 

1. k i ¥ Eb(k) is analytic on B p . 

2. There is a map k \-¥ Pb{x\ k), such that any eigenvector corresponding to Eb(k) is a multiple 
of, where Hq(k)p b (x; k) = E b (k)p b (x; k). 

3. Moreover, we can choose k i-¥ Pb(x; k), k € B p to be analytic and such that 



/ 

Jo 



\pb(x;k)\ dx 



1. 



Proof of Theorem \B.l\ L et k = fc* + k, where k will be chosen to be sufficiently small. The periodic 
eigenvalue problem (B.2I may be rewritten as 



HQ{k*)p{x; k* + k) — ( y 4nin(d x + 2-xik*) — 47t 2 k 2 ) p{x; fc* + k) = E p(x; k* + ft), 
p(x + 1; k* + ft) = p(x; k* + ft) , 



We seek an eigen-solution of (B.3l-(B.4l in the form 



p(x;k* + K) = p(x;k*) + r](x;K), 
E(k* + k) = E< + h(k), 



(B.3) 
(B.4) 

(B.5) 
(B.6) 



where we assume that r)(-; ft) _L p(-; ft). Substitution into (B.3l-(B.4| yields the following equation 
for rj(x; p, ft): 

(H(k*) - E*) Tj - (4mn(d x + 2-nik^) - 4n 2 n 2 + /j) r) = (AmK(d x + 2mk*) - 4tt 2 k 2 + p) p(-; fc«) . 

(B.7) 

Now, introduce the projection operators 

Qf = (p{-;k,),f)p(x;k,), and Q±=I-Q. 
Applying Q ± to (|B.7[) yields 



(H(k*) - E*)r] - Q± (4mK(d x + 2-KiK) - 4tt 2 k 2 + p) n = AmK Q±d x p{-; fc») = 4mKd x p(-; &*). 

(B.8) 



(B.9) 



Next, applying Q to ( |B.7[ |, i.e. taking the inner product of (B.7I with p(-; k t ), yields 
jU — 87r 2 fc,ft — 4tt 2 k, 2 + 4iri k (p(-;k»),d x r){-;p,,K)) = 0. 



We shall now solve (B.7 1 for 77, substitute the result into (B.9 1 and obtain a closed equation for 

(B.10) 



K/J.1, 77 = KTjl. 



the eigenvalue correction p = p(n). Let 

l-i 

Equations ( |B.8[ ) and ( |B.9| | become 

(H(k*) - E*)r)i - K,Q ± ('lni(d x + 2TYik i ,)-4n 2 hi + p 1 )r] 1 = 4mQ±d a .p(-;k tl ). (B.ll) 



Hi — 8n 2 k.* — 4-k 2 k + 4-Kin (p(-; d x T]i(-)) = 



(B.12) 
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Let R(E*)Q± = (H(k«) - S,) _1 Q±. Then, 
r]i(x; fii, k) = 47tz (/ — kR(E„)Qj_ (4iri(d x + 2nik+) — 4tt 2 k - 



R(E*) Qxd x p{-,k*) , 
(B.13) 



where we take |k| < p, with p chosen so that the Neumann series for the operator on the right 
hand side of (B.13 1 converges. Note that the mapping 

(pi,k) i-> rji(x; jtti,«) 

is an analytic map from {(pi, ft) : \pi\ < 1, \k\ < p'} to //p Cr (R). 
Subsitution of (B.13 1 into (B.12I gives the scalar equation 

Q{pi,n) = , 

where 

Q(pi,n) = pi — 8n 2 k f — 4ty 2 k + 4nin {p(-; fc»), d x rji(-; pi, ft)) . 



(B.14) 



(B.15) 



We now claim that (B.14 1 can be solved for pi = which is defined and analytic for 

|/t| < p , where < p < p. If this claim is valid, then r]i(x; pi(k), k) is well-defined and analytic 
in k for | ft | < p' and finally 



p(x;k,+n) = p(x;k t ) + nrji(x; pi(k), k) 
E(k, + ft) = E+ + kpi(k) 



(B.16) 
(B.17) 



are defined and analytic Floquet-Bloch eigensolutions for \k\ < p' ■ 

Now (B.14 1 is easily solved for pi = /ii(ft) via the Implicit Function Theorem. Indeed, we have 
<5(87r 2 fc„, 0) = and d lll Q(pi, K)\, air 2 k \ = 1 7^ 0. This completes the proof of Theorem 



B.l 



□ 



C The bootstrap: proof of Corollary 3.6 



We give the proof of Corollary |3.6| on the refined expansion of the bifurcation of eigenvalues of 
Hq + XV — —d 2 + Q(x) + XV(x), for Q(x) periodic. The case of Q(x) = is obtained along the 
same lines, using p (x; k ) — 1, E(k) = 4n 2 k 2 for k £ R , etc. 



Proof of Corollary 3.6 We know, by Theorem 3.4 that there exists (ip ,E ) solution of the 



eigenvalue problem (—d 2 + Q + W) ip x = E x ip x . Moreover, E x is in the gap of the continuous 
spectrum of spec(Hg) = spec(_HQ + Av), near an edge _E* = Et t (k,). In the following, we assume 
that k, — (the case where fe* = 1/2 can be treated using the same method). 

We next seek an integral equation for if) by applying the resolvent Rq(E x ) to the differential 
equation for ip x . A construction of the resolvent kernel, Rq(x,d; E x ), proceeds as follows. Intro- 
ducing the discriminant, D(E),as introduced in Appendix [A] we know that since _B* is a band 
edge and E x is in a gap, one has D(E t ) — 2 and D(E X ) > 2. Therefore, there exists ft = ft(A) > 
with 

E x = E(i\ K ) = E{-i\n), D(E X ) = e 2 ^" + e~ 2nXK > 2. 
Additionally, we define if)± = ^}±(x;E x ), the solutions of 

(-dl + Q(x))iP ± = E x i>±, ^ ± {x + l;E x ) = e ±2 * x ^ ± {x;E x ). 
More precisely, tjj± are defined as 

ip±(x) = pb t (x;^fin)e ±27vX "* x , with (CI) 

(— (d x — 2ty\k) 2 + Q(x)) pb,(x;iti) = E x p bt (x;in) , p bt (x + l;in) = pt, (x; in) . (C.2) 
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which is well-defined for A small enough, by Theorem |B.1| 

With those definitions, the resolvent operator Rq(E ) = (~d 2 + Q — E x )~ 1 has kernel 



R Q (x,y;E x ) = 



wup±] 
^+(y)ip- (g) 
w[i>±] 



if y > x, 

if y < x. 



where W[^>±] = ip' + (x)ip-(x) — tp+(x)ipL(x). Thus, for any bounded function /, 



R Q [f](x;E x 



R Q (x,y;E x )f(y) dy, 



we have (—d 2 . + Q — E x )RQ[f](x; E x ) — f. It follows that ip x satisfies the integral equation 

^ x (x) + \ f R Q {x,y;E x )V{y)i> x {y) dy = 0. 

Multiplying by Ub, (x; O)V(x) and integrating along x yields 

f V(x)u b ,(x;0)^ x {x)dx + X ff u K {x;0)V{x)R Q {x,y;E x )V(y)- l p x {y) dx dy = 0. (C.3) 
J JJm.2 

We will deduce from ( |C.3[ ) the precise behavior of k (and therefore E x — Et t (0)) as A tends to 
zero, using the following 

Lemma C.l . Let E„ — Eb„(0) be an edge of the continous spectrum, and let the hypotheses of 
Theorem 3.4 be satisfied, so that E x exists. Define Rc}(x,y; E x ) as above. Then for A > small 



enough, one has 

R Q (*,y,E X ) = Ub * ( * ; 0) e-^*-* + iff {x,y) + XkR$ (x, y), (C.4) 

47T 

where Rq is skew- symmetric: Rq (x,y) = —Rq{y 1 x); and Rq , Rq are bounded: 
\R ( °\x,y)\ + \R^(x,y)\ < Ce' 2 ^^ < C, 

where C is a constant, uniform with respect to Ak. 

In order to ease the reading, we postpone the proof of this result to the end of this section, 
and carry on with the proof of Corollary |3.6| 

By Lemma C.l and since Ub,(x;0) is uniformly bounded (see Lemma 2.51, one has the low- 
order estimate 



R Q (x,y;E x ) - 



4ty u b , (x; 0)ub, (y, 0) 



d%E(0) 2Ak 
where we used e"*"' 1 " 91 — 1 < C\n\x — y\. 



< C(1 + \x- V \+\k), 



(C.5) 



Plugging (C.5 1 into I C.3 1 and using (1 + \x\)V £ L 1 , yields 

2tt 



V(x)ub, (x;0)i(j (x)dx+ 



*d 2 k E(0) JJ R2 



u b ,(x;OyV(x)u b ,(y;0)V(y)r(y)dxdy\ < CX(1+Xk). 

(C.6) 
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3.4 



; has ||-0 A (x) — ub, (x; 0) exp (Aao|a;|) \\ LOO < A 1/4 , 
that lim J V(x)u b , (x; 0)i/j x (x)dx = J V(x)u b , (x; 0) 2 7^ 0. It follows that for A sufficiently small, 



Now we use the fact that by Theorem 



divide out J V(x)u bt (x; Q)i{? (x)dx, and deduce from ( C.6 1 
2tt 



d 2 k E(0) J R 



u b , (x; 0) V(x) dx < CXk(1 + Xk), 



from which it follows the low-order estimate of k 

I 2vr 



dim 



u b ,(x;0yV(x) dx < CX. 



(C.7) 



Let us now derive higher order estimates. For any x, y G R 2 , 

1-2-irAKla;-!/ _ ^ _j_ 2ttXk\x — y\ I < 47T 2 A 2 /t 2 jx — t/| 2 , 

so that one has from Lemma |C. II 



_Rq (x, y; E x ) - u t ,( J; 0K(i/;0)(l-2,A^- y |) _ ^ ( ^ y) 



a 2 s(o) 



Xk 



Plugging jOgll into ||C.3||, and using (1 + |a;|)V £ L 1 , yields 



< CA(l + |x| 2 + M 2 ). 

(C.8) 



Ub t {x;G)V{x)ij) (x)dx 
2tt 1 

+ 



dlE(0) K 



V(x)w6,( : r;0) 2 (l-27rAK|a;-y|)M i ,,(j/; 0)^(^)^(2/) dy 



+ ^ 7(o;)u i .(x;0)i^ 0) (a!,y)V(i/)^ A (y)dBdy| < CA 2 . (C.9) 

Let us now use that by Theorem 3.4 sup x SR \ip x (x) — u btr (x; 0) exp (Xa \x\)\ < A 1 / 4 , so 
|?/> A (:r) - u b ,(x;0)\ < C(A 1/4 + X\x\). Thus |C^]) becomes 



1 2tt 



^.(■;0)VV A ) ( 1+ „ () , 
4tt 2 

A 



K a 2 £(o) 



y v(x)«6.( 



x;0) dx 



dlE(0) 
X 



V (x)ub m (x; 0) |x — y|u(,„(y;0) V(y) dx dy 
V{x)u b ,(x;0)R { ° ) {x,y)V(y)u bt (y;0) dx dy\ <CA 1+1/4 , (CIO) 

and one deduces from |07| that J k (Ju6,(-;0)VV' A )~ 1 + 27r(a 2 £ , (0)) _1 | < CA 1/4 . Therefore, 
multiplying ^CTToj) by k (/ u 6 , (•; 0) Vi/> A ) _1 yields 



2?r 



9 2 £(0) 



V(£)m6, (x; 0) 2 dx + X 



8n J 



(dlE{0)Y 



V{x)u b , (as; 0) |as - y|uf„,(y; 0) V(y) dx dy 



X 



V(x)ubJx;0)R ( ° ) (x,y)V(y)u b ,(y;0) dx dy\ < CA 1+1/4 . (C.ll) 



Finally, we note that since Rg\x,y) — —R y Q ! (y,x) by Lemma 



(«)/ 



C.l 



the last term in (C.ll I van- 



ishes. Thus the above estimate, together with the following Lemma, completes the proof of Corol- 



lary (3.6 1. 



□ 
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Lemma C.2 . Let E* = E br (0) be an edge of the continuous spectrum, and let hypotheses of 
3.4 be satisfied, so that E x exists. Then for X small enough, one has E x = E(iXk), and 



Theorem 
E x — E, 



,X 2 K 2 d 2 k E b ,(Q) + 0(X A ) 



Proof. We Taylor expand D[E) about E„ = E b , (0). 

D{E) = D(E t ) + D'{E*){E -E*) + ((E - E,) 2 ) = 2 + D'(E„)(E -E*)+<D ((E - E,) 2 ) . 

(C.12) 

Let's first apply ( |cl2| to E = E K (fc) in the spectral band. One has D{E b , (fc)) = e 27,ik + e~ 2vik = 
2 -4n 2 k 2 + (fc 3 j~Fnially, since d k E b ,(0) = d A k E b , (0) = 0, one has E b ,(k) = E, + ±<9 2 £(0)fc 2 + 
O (fc 4 ) . Identifying with ( [cl^ ), it follows D' {E,){\djE{Q)) = -4tt 2 , thus D'(E,) = g ^f W ■ 
Next let's apply jCHg) to £ = £\ recalling D(E X ) = e 2 ^« +e - 2 ^ = 2+4vr 2 A 2 ft 2 +O (A 4 ft 4 ) . 
One has from p3| in Theorem that E x - E« = 0{X 2 ), and from ( [07) that /s = 0(1). 
Consequently, ( C.121 yields 



, 2,22 

47T A ft 



D'{E«)(E X - E,) + 0(X 4 



d 2 k E b ,(0) 



(E x -E b ,(0)) + C(A 4 ). 



Finally, we deduce E x - E* = -~\ 2 K 2 d%E bt , (0) + C(A 4 ) and the lemma is proved. □ 

We conclude this section by the proof of Lemma |C.1| 
Proof of Lemma \C.l\ Let us Taylor-expand ip±, as defined by ( |C.1[ )- |C~2"| ). By Theorem |B.1| one 

has ip±(x)e ll2lTXKX = p bzt (x; =FiAft), thus 



i'+(x)e = p b ,(x-0) - iXndkPb, (x; 0) 

7/)_(x)e 2 ' rA ' la: = p&„(a;;0) + i\nd k p bt (x; 0) 



^a 2 p(x;0) + i^fyt), (C.13) 
%^a 2 p(s;0) - i^!^p(a;;< 7 _), (C.14) 



with —Aft < 7+ < < 7_ < Aft. 

Remark C.3. ./Vote that by the construction of p b (x;k* + ft) in Theorem B.l we have that k 



p b (x;k* + ft) £ L (R) is analytic in a cornple 



•orhood ftj < fti. By t/ie equation for p b , 



ft h-> p b {x;k, + ft) £ i/ 2 (R) is analytic and thus d k p b (x;k) and d x d k p b (x;k) are well-defined and 
uniformly bounded for k near fc* and x in any compact set. 

Since p bt (x; 0) = u bt (x; 0), it follows 



W&±]R Q (x,y;E X ) = 



(u b , (x; 0)u K (y; 0) + iAftr (0) (x, y; Aft) + (Aftfr^ (x, y)) e 27TXK{x ~ v) if y > x, 
(u b ,{y;0)u b ,(x;0) + iXKr (o) (y,x;XK) + (Aft) 2 / 1 ^, y)\ e 2 ^(»-*) if y < a:. 



(C.15) 



with 



r (0) (ir,y; Aft) = p b , (x; 0)d k p br (y; 0) - d k p bst (x; 0)p bt (y; 0) = -r w (y, a; Aft), 



(«)/ 



and r±'(a;,y) is bounded, uniformly with respect to Aft. 

Let us now turn to VK[^±] = tp' + (x)ip-(x) — ip+(x)il>'_(x). From ( C.13 H C.14 1, one has 

W[i>±] = 2Aft( 2ivp br {x;0) 2 - ip bt (x; 0)d x d kPb , (x; 0) + i{d x p b , (as; 0)) {d kPb , (x; 0)) \ 

+ (Xn) i w r (x; Aft), 

with u)r(ic) uniformly bounded, independently of x and Aft. 
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Since W^i] is independent of x, one has 

W[i>±] = f W{i>±] 
Jo 

= 2Ak J ( 2np bt (x; 0) 2 - ip b , (x; 0)d x d k pb, (x; 0) + i(d x Pb, (x; 0)) (<9 fc p&, (x; 0)) ) da: 

3 f 1 

+ (Ak) / TO r (a;;Afi;) dx. 
Jo 

Using that pb, [x; 0) 2 da; = u b , (x; 0) 2 da; = 1, one deduces 

W[ip±] = 2\k(2-k + 2i f p(x;0)d x d k p{x;0) dxj + O ((Ak) 3 ) . (C.16) 



Now, let us recall that p b , (x; in) satisfies ( C.2 1 . Deriving twice with respect to k — in, one obtains 



(—(d x — 2ttk) 2 + Q(x) — E{in)) dlp(x; in) = 2dkE(in)dkp(x; in) + dlE(in)p(x; in) 

— 8ni(d x — 2nn)dkp(x; in) ~ 87r 2 p(a;; in) . 

We now apply this identity at k = 0, and take the inner product with p{,„(a:;0). It follows 
= 9fc.E(0) - 8iri p(x; 0)d x d k p(x; 0) da; - 8tt 2 . Therefore, {CJL61 becomes 



WW>±] = 2Ak ^ jE(0) + 0((Ak) 3 ). (C.17) 

47T 



Finally, ( C.15 1 and ( C.17 1 clearly imply (C.4|, and Lemma |C.1| is proved. □ 
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